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A closed-form formula is derived for the generalized Clebsch-Gordan integral /3j [P v (x)] 2 P v (-x)dx, with P v being the 
Legendre function of arbitrary complex degree v e C. The finite Hilbert transforms of P v (x)P v (-x),-l < x < 1 and 
(1 + x) v ~ n P v (x), -1 < x < 1 are evaluated for any non-negative integer n satisfying Re(v-n) > -1. An analytic proof is 
provided for a recently conjectured identity Jq[K(V 1 - k 2 )] 3 dk = 6 / 1 [K(fe)] 2 K(v / l -k 2 )kdk = [r(±)] 8 /(1287r 2 ) involving 
complete elliptic integrals of the first kind K(&) and Euler's gamma function T(z). 

Introduction 



In a wide range of mathematical and physical contexts, one needs to evaluate multiple elliptic integrals where 
' the integrands involve one or more factors of complete elliptic integrals. As a glimpse of some recent applications in 
number theory, high-energy physics, statistical mechanics and probability theory, we mention automorphic Green's 
functions [1], multi-loop Feynman diagrams lattice Green's functions [3] and short uniform random walks 
(J. Various techniques, of algebraic, geometric, combinatoric or analytic flavor, have been developed to express these 
multiple elliptic integrals in terms of familiar mathematical constants and special values of Euler's gamma function 
fli2,3,4,5,i. 

(T) | When the integrands contain the products of three or more complete elliptic integrals, the evaluation can become 
■ analytically challenging. For example, the following relation 



[r(i)] s 



1287T 2 



= J" 1 [K(v/l-/e 2 )] 3 dA! = 6 pmkrfK^Vl-k^kdk 



has been recently conjectured by J. G. Wan in [6] based on numerical evidence. However, a rigorous proof seems to be 
■ beyond the available toolbox for elliptic integrals. 

In this work, we solve Wan's conjecture and compute a related family of multiple elliptic integrals with sophisticated 
T 7"! ■ appearance. We employ an approach that is probably overlooked in previous literature: connecting elliptic integrals 
J> ' to spherical functions and spherical rotations. Concretely speaking, the complete elliptic integral of the first kind 
sin 2 9) 1/2 dd is related to fractional degree Legendre functions P_i/2,P-i/3,P-i/4,P-i/6, where 



, (2v+l)^ 



c3 • 2 r° cos- , 

P v (l) = l; P v (cos0):=- g dp, Be(0,n),v£C 

xJa ^(cosjS-cose) 

extend the familiar Legendre polynomials 

1 d^ 

p ' (x)= 577fT-? [(r2 - 1) ' 1 ' £eZ ^° 
2 e £\ dx l 

to arbitrary complex degree v e C. By a modest generalization of the "spherical harmonic coupling" for Legendre poly- 
nomials (Clebsch-Gordan theory) to Legendre functions of arbitrary degree, we are able to prove the integral identity 



\ pni2p , v l + 2cos(;r Z ) nT(^)T(^) 
[P v (x)rP v (-x)dx = hm - , veC 

3 [r(^)f[r(f)fr(f) 



and compute the finite Hilbert transform 



^ r 1 2P v (()P v (-0 , , [Pv(x)] 2 -[P v (-x)] 2 

& I ; « — df = — — , veC\Z. 

J-l 7i(x-i) Sin(V7T) 
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The pair of formulae above, along with a couple of geometric transformations related to rotations on spheres, allow us 
to evaluate a slew of multiple elliptic integrals, which embody Wan's conjectural identity as a special case. Perhaps an 
intriguing twist, the finite Hilbert transform of 2P V (£ )P v (-<f) is also equivalent to the following identity: 

9 f 1 (1 + O v ~ n Pv(OT^-r = (1 + x) v - n lim - \cos(nz)P z {x) - P z (-x)l, n e Z> ,Re(v -n) > -1. 

J-l 2(x-$) z—*y 2sin(7rz) 

The article is organized as follows. In |TJ we review some classical results that express fractional degree Legendre 
functions as complete elliptic integrals, and give a new proof of Hobson's coupling formula for two Legendre functions, 
based on scaling limits. We carry on the scaling analysis in so as to interpolate the standard Clebsch-Gordan in- 
tegrals for Legendre polynomials into a closed-form expression for f\ [P v (x)] 2 P v (-x)dx with arbitrary complex degree 
v £ C. This is followed by a geometric intermezzo on spherical rotations in 5j3l providing tools to transform one multiple 
elliptic integral into another, as well as preparing some integral identities for 5j4] The finite Hilbert transforms (also 
known as the Tricomi transform) of the functions P v (x)P v (-x) and (1 + xY~ n P v (x),n e Z>o,Re(v-re) > -1 are studied in 
|4j in the spirit of "spherical harmonic coupling". Lastly, serves as a perspective on subsequent works on multiple 
elliptic integrals and their applications. 



1 Fractional Degree Legendre Functions and Hansen-Heine Scaling 
Limits 



The Legendre functions of the first kind can be defined via the Mehler-Dirichlet integral 

(2v+l)/3 
5 2 

v/2(cosj8-cos0) 

Alternatively, we may characterize P v (x),-1 < x < 1 as the unique C 2 (-l,l) solution to 

d 



2 ,e cos^HM 

P V (D=1; P v (cos0):=-| 2 — dp 0e(O,7r),veC. (1) 

nJo \/2(cos/3-cos0) 



dx 



(1-x 2 )^ 
dx 



+ v(v +!)/•(*) = 0, /XI) =1. (2) 



By either definition, we have P v (x) = P_ v _i(x). Later in this work, we shall also use the Legendre functions of the 
second kind, given by 

71 

Qv(x):=—— — -[cos(vjr)P v (x)-P v (-x)], veC\Z; Q n (x) := limQ v (x), neZ> (3) 
2sm(v7r) v—n 

for -1 < x < 1. By convention, the Legendre function of the second kind Q n (x) is undefined for any negative integer 
degree n e Z<o. While both P v (x) and Q v (x) solve the same Legendre differential equation of degree v e C \ Z<o, the 
function Q v (x) blows up logarithmically as x — » 1 - + , unlike the natural boundary condition for P v (l) = 1. 

Additionally, the Legendre function of the first kind can be represented as the Laplace integral under some restric- 
tions: 

P v (z):=— e vheU+iVTZ ^ C0S(l>) dcl>= — eViogfe-ivCTco80) d ^ R ez >o,veC. (4) 
2n Jo 2n Jo 



1.1 Kleiber-Ramanujan Correspondence 

Highlights of the relation P_i/2(cos0) = (2/;r)K(sin(0/2)) date back to J. Kleiber's posthumous publication in 1893, 
where systematic studies were carried out on connections between complete elliptic integrals and Legendre functions 
P(2„+i)/2,« e Z, the latter of which are useful for solving differential equations in toroidal coordinates [Ref. 8, §§253-258]. 

In his seminal paper ||9J] and legendary notebooks [Ref. 10, Chap. 33], S. Ramanujan developed rich and elegant 
theories relating P-i/3,P-i/4,P-i/6 to complete elliptic integrals of the first kind. 

We recapitulate the aforementioned classical results in the following short lemma. 

Lemma 1.1 (Some Fractional Degree Legendre Functions) (a) We have the following identities for < 9 < n: 



2 [ B\ 2 r M/2 dw 

rr i _ _ _ i — , 

2 > nJo Vl-sin^e^sin 2 !//' 



P_ 1/2 (cos 0)= -Rising] = - P ==. (5) 

Ji { 2j n Jo 



„ , „ s 2 1 I 2sin(0/2) 1 2 r M/2 dw 

P_i/ 4 (cos0)= K \ . \ n ' \ = - I y (6) 

n Vl + sin^) iV l + sin(0/2)j n Jo ^l-si^^sinV 
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(b) For a parameter p e [0, 1), the following relations 
21p 2 (l + p) 2 \ /l + p+p 2 . 



-1/6 



-1/6 



1-2 



4(1 + p+p 2 ) 3 
27p 2 (l + p) 2 



4(1 + p + p 2 ) 



2l3 



1 = 



l + 2p 

1 + p+p 2 
l + 2p 



-1/2 



-1/2 



ll-2 P(2+P) ) 


= l\ 


jl + p+p 2 K 




\ l + 2p J 




1 l + 2p 




( 2 f + ^-l] 


= l\ 


jl + p+p 2 K 




I l + 2p J 




1 l + 2p 





lp(2 + p)\ 



l + 2p j 



1-P 



2\ 



l + 2p) 



hold for degree v — - 1/6, and the transformations 



-r-l/3 
P-V3 

apply to degree v 



1-2 



27p 2 (l + p) 2 
4(l + p + p 2 ) 3 

, 27p 2 (l + p) 2 x 
' 4(1 + p+ p 2 ) 3 

-1/3. 



1 + p+p 2 

VTT2p~ 
l+p+p 2 



V3 + 6p 



-1/2 



-1/2 



1-2- 



p 3 (2 + p)^ 2 1 + p + p z 
l + 2p 



, p ,j (2 + p) 
l + 2p 



n v/l + 2p 

2 1 + p + p 2 
^ ^3 + 6p 



K 



p 3 (2 + p)\ 



l + 2p j 



K 



(l-p)(l + p) 
l + 2p 



3\ 



(7a) 
(7b) 

(8a) 
(8b) 



Proof (a) Let K(k) = Jq' 2 Vl-k 2 sin (f>d(p be the complete elliptic integral of the second kind, then we have the differ- 
ential relations 

dK(A0 E(k) K(k) d~E(k) E(k)-K(k) 



dk k(l-k 2 ) 



dk 



It is thus routine to check that the proposed forms of elliptic integrals in Eqs. [5] and [6] satisfy the said Legendre 
differential equations and natural boundary conditions (Eq.©. 

The last integral in Eq. [5] is the standard definition for the complete elliptic integral of the first kind. With a 
variable substitution Q = arctan( y/1 + sin(0/2) tan y/) [Refill], p. 105], one can convert the last expression in Eq. |6]into 
the standard form. 

(b) As in the proof of part |(a)[ one can verify the relations between the respective fractional degree Legendre functions 
and complete elliptic integrals by a routine, though time-consuming procedure of Legendre differential equations. 

Equations[7a]and[7b]have appeared in Ramanujan's work [Ref. 10, pp. 161-163] on the elliptic function theory for sig- 
nature 6. As suggested by [ 10], Eqs.[7al[7b]can be constructed from first principles, using a hypergeometric transforma- 



tion that identifies P_i/6(cos0) = 2F1 [ 



1/6,5/6 1 
1 



1/12,5/12 1 



sin 2 0),O: 



I sin 2 |) ,0 < 8 < nl2 with the Fricke-Klein function 2-F1 ( , 
9 < n/2 [Ref. [13, P- 334], the latter of which has a well-known connection to complete elliptic integrals of the first 
kind H. 

We may combine Eqs. [7aland[7bl into a single formula 



-1/6 



x(9-x 2 ) 



(3+x z y 



2\3/2 



v^Tx 2 2 
— - — P-m(x) = - 
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v / 3 + x 5 



K 



ate (-1,1], 



(9) 



a form that is handier for conversion between integrals involving P-i/e and K. 

Ramanujan's studies of the elliptic function theory for signature 3 bring us Eqs. [8a] and Hb] One may consult [Ref. 
pp. 112-114] for the detailed theoretical arguments leading to the discovery of the algebraic relations between P_i/3 
and K. ■ 

According to the recursion relation for Legendre functions (x 2 -l)dP v (x)/dx = vxP v (x)- vP v -\(x) and the symmetry 
P v (x) = P- v -iix), one can relate some fractional degree Legendre functions to complete elliptic integrals of the first and 
the second kinds, such as 



P 1/2 (cos#) = 



2E sin- 



K sin- 



P 3/2 (cos0) = 



Pl/4(COS0) = — 
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2v / l + sin(0/2)E 



2sin(0/2) 



1 + sin(0/2) J 



3tt 
1 



8cos0E 
I 



( ■ 6 ) 
sin — 

\ 2) 



(4cos0 + l)K 



( ■ d ) 
sin - 

\ 2) 



P 3/4 (cos0) = — 
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2v / l + sin(0/2)E 



2sin(0/2) 1 



1 + sin(0/2) 



Vl + sin(0/2) 
l-2cos0 



2sin(0/2) 



1 + sin(0/2) J 



Vl + sin(0/2) 



2sin(0/2) 



1 + sin(0/2) 
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1.2 Hansen-Heine Scaling Limits and Hobson Coupling Formula 

By taking scaling limits of certain identities for Legendre polynomials, one may extend their validity to Legendre 
functions of arbitrary complex degree. In this manner, one can compute some multiple elliptic integrals by interpolating 
Legendre polynomial identities, the latter of which often have clear motivations from the harmonic analysis on spheres. 

As we may recall, for positive integers m, the associated Legendre functions of degree v and order m are defined by 



p;(x):,(-ir(i-, 2 r /2 ^, 

ax m 

Q v (z).-(z -1) , 



T(v + m+l) 



One also writes P° = P v and Q® = Q v . The spherical harmonics 

Y em (6,<p):= 



2£+l(£-m)\ 



An (£ + m)l 



Pf{cos8)e imip , CeZ> ;mEZnKfl 



carry the Condon-Shortley phase Y (m (0,<p) = (-l) m Y e _ m (0,0). 

In the following, we show that the studies of Legendre functions of arbitrary degree can be built on the Mehler- 
Dirichlet formula for spherical harmonics (see item 8.927 in 11311 or §4.5.4 in 1 14]), which couples two normal vectors 
»i(0i,0i)= (sin 0i cos sin 0i sin 0i, cos 0i ) and «2(#2,02) = (sin02COS02,sin02 sin</>2,cos02) on the unit sphere: 



- -£P/(ni-n a )cofl[(/+|)A|: 



4»I E 

e=0m=-e 



2(+l 



0ij(cos/3-rai -ra2) 
^(cos/S-ni -712)' 



|»ll = |n 2 | = l. (10) 



Here, the Heaviside theta function is defined as 8h(x) = (1 + ifr)/2,je sR\ {0}, and the range of the angle p is (0,;r). 

Lemma 1.2 (Hobson Coupling Formula) For 61,82 e [0,7r),v e C, i/ie Legendre function of the first kind satisfies the 
Hobson coupling formula 



1 f 2jI „ 
2k Jo 



(cos 0i cos 02 + sin 0i sin 02 cos <p)d<p- 



P V (COS01)P V (COS0 2 ), 01 + 2 <7T 

P v (-COS0i)P v (-COS0 2 ), 01 + 2 >7T 



(11) 



which specializes to the following integral identities for 61,82 G [0,7r): 

'K(sin^)K(sinf ), 0i + 2 <7r 
k(cos^)k(cos%), 0i + 2 >tt 



1 C 2n 



-A 

4 Jo 



K[sin^]d0 = 



2sin(0/2) 



+ sin(0/2) 



d0 



1 jt( / 2sin(6i/2) 

^l+sin(0i/2Vl+sin(0 2 /2) IV l+sin(0i/2) 



v 7 ! + sin(0/2) 



2cos(0i/2) 



Vl+cos(e 1 /2)Vl+cos(0 2 /2) IV l+cos(0i/2) 

where cos& = cos0i cos02 + sin0i sin02 cos0. In particular, we have 



(12) 



(13) 



K|sin- 



K sin- K cos- = 



l + \/u 



K 



( / 2A 



l + \/u 



V2 
1 + Vu 



K 



f / 2y^ \ 
1+Vu 



K 



1 + y/u 



run 

\ K(sin0cos0)d0, 
Jo 

J K^l-sin^cos^jd^, 
j^-gi I 4Vu(l-ujcos(p 

K 

Jo 



We[o,|], C2f) 
V0£[o,|], fi2f) 



1 + 2\/u(l- u)cos0 



d0 



/ 2v / 1-4m(1-m)cos 2 1 



V / l + 2 v / w(l- u)cos< 
d0 



1 

°'2 



1+ \/l-4u(l-u)cos 2 y 



Vl + v / 1-4m(1-m)cos 2 



G3f) 
Vwe(0,l). CEHf) 
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Proof When v = n is a non-negative integer, one can verify Eq.[Tl]by resorting to Eg. 1101 



1 f 2jI 

— | P„(cos0icos02 + sin0isin02costf>)drf> 
2n Jo 

f2n \ 2 ft ff (cosj8-cos#icos02-sin0isin02COS(/>)cos[(7i+ |)/3]d/3 



l r Z7T 2 r 
2tc Jo ti Jo 

2n Jo n Jo e = o m =-e 



yj cos p - cos 0i cos 62 - sin 0i sin 62 cos </> 
' Y> n ,(0i,O)Y /m (e2,0)cos[(/+ |)j8]cos[(«+|)j8] 



d0 



2^ + 1 



d/3 



d0 = P„(cos0i)P n (cos0 2 ), 



where spherical harmonic modes other than Y n o will have vanishing contribution upon integration in the angular 
variables (5 and <p. Here, due to the parity relation P n (-x) = (-l) n P n (x) for non-negative integers n, one always has 
P n (cos0i)P n (cos02) = P n (-cos0i)P n (-cos02), irrespective of the sum 61 + 62- 

Fixing two positive parameters s and t, one may consider a neighborhood of the origin in the complex z-plane such 
that Recos(sz) > 0,Recos(£z) > 0,Re[cos(sz)cos(iz) + sin(sz)sin(fz)cos0] > 0. Then, given an arbitrary v £ C, we can 
apply the Laplace integral representation for the Legendre functions (Eq. [4} to evaluate the following scaling limit: 

r 2n dch 

fs t-v(z) := P v / z (cos(sz))P v / z (cos(£2:))- / P v / (cos(s2:)cos(£2:) + sin(sz)sin(£2:)cos0) — 

Jo 2tc 

-f 

Jo 



c2tt 



■ log(cos(sz)+z sin(sz)cos(fti ) *^01 f ^ log(cos(£z)+z sin(tz)cos0?,) ^02 

2n Jo 2ji 



f \ f log(cos(sg)cos(f2)+sin(s2)sin(^z)cos0+?\/l-[cos(s2)cos(^2)+sin{g2)sin(^)cos0] 2 cos0^) ^03 1 d0 

Jo iJo 2it J 2n 



r2n 

Jo 



e 

? ivscos(/)i ^01 

2n 



r-2n 

Jo ' 



d02 
2n 



r-2n i n2n 

Jo I Jo 



V / s 2 + f 2 -2sico S </.co S( />3 1 £0 



2n 2n 



0, as z^O, 



where the last equality turns out to be a special case of the Sonine-Gegenbauer formula for Bessel functions [Ref. I IE 
§12.1]: 



1 C 2n i 

Jo(vs)Jo(vt)= — \ Jo(vV s 2 + t 2 -2sfcos0)d0. 
2n Jo 



Thus, with fixed parameters s, t and v, the function f s ,t-v(z) is analytic near the origin in the complex z-plane, and 
fs,t;v(z) = on the set {0} U {vln\n £ Z n [0, +00)}, which contains an accumulation point {0} in the domain of analyticity 
for f s ,t;v(z)- Therefore, the function f s ,t;v(z) must vanish identically in a certain neighborhood of the origin, the spatial 
extent of which depends on the nature of the positive parameters s and t. 

Now suppose that s + t <n, then cosscosi + sins sin t cos (p > cos(s + t) > -1. Accordingly, for a complex number z 
sitting in a certain open neighborhood of the unit interval [0, 1], the expression cos(sz)cos(£z) + sin(sz)sin(iz)cos0 will 
miss the value -1, which is the logarithmic branch point of the Legendre functions P^ of non-integer degree fi £ Z. 
Hence, when s > 0, t > and s + t <n, we have a univalent complex-analytic function f S! t;v(z) in an open neighborhood 
of the unit interval [0,1], where resides a non-isolated set of zeros. By the principle of analytic continuation, we may 
conclude that f s ,t,v(X) = for s > 0, t > and s + t <n, which extends, by continuity, to the s + t - n scenario as well 
as those situations where st = 0. This proves Eq. [11] for the cases of 61 + 62 £ n. The rest of Eg. [TT1 follows from the 
invariance of the integral 



2nJo 



P v (cos 61 cos 62 + sin 0i sin 62 cos 0) d 



under the equatorial reflections (61,62) —* (n - 6\, n - 62). 

Setting v = -1/2 in Eq.Qj] we o btain Eq.[12l setting v = - 
Eq.[12](resp.[T3]! becomes Eq.[l|l](resp.[llt]l. When 6 1 = n- 
procedure on Eq. [13] would bring us the evaluation 



1/4 or v = -3/4 in Eq.[TT] one verifies Eq. [T3"| Wh en 6 l = 6 2 , 
-62-6, one may specialize Eo. 1 121 into Eq. ll2n A similar 



r n/2 

K 

Jo 



2v / 1-4u(1-m)cos 2 1 



1+ v / l-4u(l-u)cos 2 



V / l + v / l-4u(l-u)cos 2 



Vi + VuVi + VT- 



2v/u \ 



K 



2VT- 



l + s/T- 



V M e|0,- 
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while the substitution u = (1 - t) 2 /(l + t) 2 and Landen's transformation [Ref.[H, item 8.126] allow us to see that 

2Vt 



2 V / T :r w 



j l + Vl-u l 



l + y/t~[l + Vt 



K 



l + y/u 



K 



1- Vu 



1+ Vu, 



This verifies Eg. ll3n for u e (0,1/2]. Furthermore, the right-hand side of Eg. [13T l is intact under the transformation 
u« 1 - u, whereas the symmetric extension of its left-hand side hinges on the identity 



l + Vu 



1- Vu 
l+Vu 



K 



1 + \fu 



i + 



l-VT- 
l + VT- 



K 



2^/1- 



1+ VT-Uj 



0<u< 1, 



which is an equivalent formulation for the products of two Landen's transformations: 



1 + t 



i+t 



2 " "' (1 + Vt) 2 \1 + Vt 

with the correspondence of variables being u = (1 - t) 2 /(l + 1) 2 . 



K 



1 + v 7 * 



0< t < 1, 



(14) 



Remark (1) Using Eq. [TO] we may readily adapt our proof of the product formula in Eq. QT]to associated Legendre 
functions: 



1 r(v + m + l) f 2n „ ,„„.„.„ 

/ P v (cos0icos02 + sin0isin02COS(/))cos/n</)d0 = 

: + 1) Jo 



2n T(v-m- 



P™(cos0i)P™(cos<9 2 ), e x + e 2 <n 
P™(-cos0i)P^(-cos0 2 ), 01 + 02 >K 



en?) 



where 0i,#2 e [0,tt), v £ C, m £ Z. Here, the corresponding scaling limit is the Sonine-Gegenbauer formula 

f27T 



l r Z7t I 

J m (vs)J m (vt) - — I Jo(vvs 2 +t 2 - 2st cos (b) cos md) deb. 
2n Jo 



(Doubtlessly, the cases of 



P^ 1/2 (cos0): 



;rsint 



E I sin - I - K I sin - I cos z - 



Pl 1/4 (cos0) = P^ 3/4 (cos0): 



v / T+sin(072) 



E 


/ 2sin(0/2) ' 


-K 


' / 2sin(0/2) ' 






V 1 + sin(0/2) j 






V l + sin(0/2) 







will then lead to generalizations of Eqs. [12] and [13] into integral formulae involving complete elliptic integrals of the 
second kind.) Effectively, we can combine Eq. 1111"! with Eq. [11] to develop a Fourier expansion in (f> for the function 
P v (cos0icos#2 + sin0isin#2cos0) (necessarily under the constraints < 9\ < n,0 < 62 ^ Ji,0 < 61 + 82 £ n), which turns 
up as a uniformly convergent series formerly known to E. W. Hobson [Ref. 8, §226]: 

°2 T(v-m+l) 

P v (cos0i cos02 + sin0i sin0 2 cosrf>) = P v (cos0i)Py(cos0 2 ) + 2 Y P™(cos0i)P™(cos0 2 )cosm0. 

Hence we refer to Eq. Qj] as the Hobson coupling formula. Hobson's original approach draws on the Fourier series of 
(z + Vz 2 - lcos</>) v for complex-valued v, which is not a method based on scaling limits. 

(2) The representation of Bessel functions J m as the scaling limit of (associate) Legendre functions P™ is a special 
case of Hansen's formula rRef. flH. §5.7]. Heine extended the limit procedure so that one may connect Bessel functions 
Y m and the (associate) Legendre functions in a similar fashion [Ref. 15, §5.71]. The Hansen-Heine scaling limit 
procedure (m = 0) will be used again during the proof of Eq. H9j„ „)1 namely, 



[P v (x)] 2 P v (-x)dx^hm 1 + 2COS(7r2) 

2— v 3 



»r(f)r(^) 



[r(^)] 2 [r(^)] 3 r(^)' 



veC 



in Proposition [2TT] 

(3) The integral formula in Eq. ll2fl is usua lly attributed to Glasser [16] , and has been mentioned recently in the work 
of Bailey et al. J2]. We will recover Eq. ll2fl later in Proposition l3.3l (see S3. 2D . using a method independent of the Hobson 
coupling formula for Legendre functions. 
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(4) The Hobson coupling formula is sensitive to the geometric constraint < #i + 02 ^ n. Thus, it is not possible to 
directly extrapolate the formulae in this lemma to an integral representation for [P v (x)] 2 , -1 < x < where v is arbitrary. 
Practically, this obstacle can be overcome in various ways. For v = -1/2,-1/3,-1/4,-1/6, the angular restriction in the 
Hobson coupling formula can be circumvented by using Ramanujan's integral representation for [K(/e)] 2 ,0 < k < 1 (see 
Proposition 13.31 1. For generic v, Hobson's approach yields an integral representation for P v (x)P v (-x),-l < x < 1 and 
[P v (x)] 2 ,0 < x < 1, and the finite Hilbert transform in Proposition 14.21 (see 34.21 1 will enable us to express [P v (x)] 2 - 
[P v (-x)] 2 ,0 < x < 1 in terms of P v (<f)Pv(- 0,-1 < {< 1. □ 



2 Generalized Clebsch-Gordan Integrals 

As the C 2 (-l, 1) solutions to the Legendre differential equation of arbitrary degree v e C \ Z<o: 

+ v(v +l)f(x) = 



_d_ 

dx 



a-x 2 )^ 

ax 



(15) 



are exhausted by linear combinations of P v (x) and Q v (x), a special case in Appell's theory of third order ordinary 
differential equations 111711 then tells us that the C 3 (-l, 1) solutions to 



_d_ 

dx 



(1-x 2 ) 



dx 



(1-x 2 )^ 
dx 



+ 4v(v + 1)(1 - x')f(x) \ + 4v(v + l)xf(x) = 0, v £ C \ Z< 



(16) 



belong to a three-dimensional space spanned by [P v (x)] 2 , P v (x)Q v (x) and [Q v (x)] 2 (or by an equivalent basis set 
{[P v (x)] 2 ,P v (x)P v (-x),[P v (-x)] 2 } when v t Z). 

The "interpolation" from Legendre polynomials Pe(x), £ e Z>o to non-integer degree Legendre functions P v (x), v £ Z 
and the "scaling limit" for large degrees will enable us to extend some Clebsch-Gordan integral formulae to generic 
Legendre functions, as elaborated in the proposition below. 



Proposition 2.1 (Generalized Clebsch-Gordan Integrals) Define 

Tfi, v ~ j ^P f i(x)P v (x)P v (-x)dx, M,veC, 

then we have the symmetry 



the recursion relation 



0. + + I) 2 - (2v + D 2 ]T M+1 , V - A, 2 ~ (2v + ifW.-i,. = 4(2 M+ l)sin(^)sin(v,) ; 



(17) 



(18) 



the representation in terms of generalized hypergeometric series 



2 sin(^7r)sin(v7r) 

7I 2 



-4F3 

V 



I 1 IzE Ell _ v 
-*•> 2 ' 2 ' 



V 2 ' 2 ' 



1-V 



V+ 1 



1^3 



2 ' 2 ' 

Hzif M±2 v + 2 
2 , 2 > v + z 



2 sin(/ijr)cos(v7r) 



— 5-^4 



/ 1 I-E-v i + v 

2' 2' 2' W '- L ^ v 



1 2-£ l-2v 2v+3 
V 2 ' 2 ' 2 



/i + 1 



/ll 1+u 1 , 

' 2>2> 2 ,-V,l + V 

1 3+£ l-2v 3+2v 
-"-> 2 ' 2 ' 2 



(19) 
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along with some generalized Clebsch-Gordan integral formulae 

nT [ 2 ^y[ 2 ^) 



r (^)] l r (fl] r(^)r(^) 
71 cos(v^)r ( 2v+1 2 - 2m ) r(v + 1 + m) 



, fiCZ,neZ>o 



m e Z>o,v+ • 



[r (i=fg)] 2 [r(m+i)] 2 rev +1 - m)r ( 2y+ | +2 '" ) ' - - 2 



(-1)" [T(m + ^)] 2 r(m -n - \)I(m +n + §) 
^ [r(m+l)] 2 r(m-re)r(m + « + 2) ' 



to e Z>o,«- EZ,m — n€ Z<o,to + n + 2 £ Z<o 
to £ Z>o,(to - re e Z<o or to + 7i + 2 e Z<o) 



flSU»)) 

G9j2 m , n+ i)) 



7 



T vv - lim 

' z—v 



i + 2cosfa z ) gr(^)r(3f2) 

3 [nW^mfn^r 



veC 



_ ,. sin(7rz)sin(27rz) 
T 2 v-i,v=lim 5-5 , veC 

_ ,. sin(n z)sin(2 jiz) 
T2v+2,v = ~ lim 7T. 7-5 , veC. 



GUv.v)) 

dfev-l.v)) 
G§l2v+2,v)) 



Proof The symmetry of Legendre functions with respect to degree P v (x) = P_ v _i(x),Vv e C explains the related prop- 
erty of T^ v (Eq. [1711- 
Judging from the familiar recursion relations for Legendre functions 

and the differential equations given in Eqs.[15]and[16l we may carry out the following computations: 
v(v + + l)T M+ i iV + ^T M _i, v ] 
= v(v+l) J [(^+l)P f i + i(x) + ^P fl -i(x)]P v (x)P v (-x)dx = (2p+l)v(v+l) J xP M (x)P v (x)P v (-x)dx 

- 2 " +1 fV>.± 

J-i dx 



4 

4 J-i 
2p+l f 1 
4 

2/*+l 



(1-x*)— 
dx 



2 d(P v (x)P v (-x)) 
(1-x )- 



dx 



a-x A ) 



dx 



„ 2 d(P v (x)P v (-x)) 
(1-x )- 



(1-x 2 )— 
dx 



dx 

„ 2 d(P v (x)P v (-x)) 



+ 4v(v + 1)(1 - x^)P v (x)P v (-x) } dx 
dP fl (x) 



dx 



+ 4v(v + 1)(1 - x z )P v (x)P v (-x) 
dP M (x) 



dx 



■dx 



dx 



- dx + n(n + l)v(v + l)(T„_i v - T„+i v ) 



M(jU+l) (l-x z )P A ,(x) 



d(P v (x)P v (-x)) 4sin(/i;r)sin(vjr) 



dx 



-dx 



+ + Dv(v + lXT„_i, v - r„+i, v ), 



where we have used in the last line the properties P v (l) = 1, Vv e C and 



2 ,dP v (x) 2 sin( vn) 

hm (1-x) = , VveC. 

-1+0+ dx n 



We may then proceed with the simplification 



(2ft+l) 



£ (i -* 2 



)P M (x) 



d(Py(x)Py(-X)) 

dx 



dx 



/1 <jp ( x ) />i 

(l-x 2 )Py(x)P v (-x) — r — dx + 2(2u+l) / xP„(x)P v (x)P v (-x)dx 
1 dx J-i 

= lj(h + 1XT m+ i,v - Vi,v) + 20* + DT M+ i,v + 2ju7y i )V) 



so as to confirm the recursion formula of T M>V (Eq. 

For non-negative integers m and n meeting the requirement m < 2n, we might recall from the standard Clebsch- 
Gordan theory that 



p m (x)p n (x)p n (-x)dx=2(-irr n n \ =(-ir — - — 2 *■ 2 2 9 2 1 9 

-i y° °) r ^ 2 r 2 r r ^ 



QSUi.)) 



where ^ 2 2 ^ 3 3 ) is the Wigner 3-j symbol. The right-hand side expression in Eq. \19\„ „M s regarded as zero if m 
is a positive odd integer, as anticipated from the relations P m (x) = -P m (-x) and Ti^ 1 ) = oo in such circumstances. 
We have T m ,„ = f^ 1 P m (x)P n (x)P n (-x)dx = for non-negative integers m and n satisfying m > 2n, as indicated by the 
respective poles of the gamma factors in the denominator on the right-hand side of Eq. [T9l m „\\ 

Actually, even without the prior knowledge of the standard Clebsch-Gordan theory, one can start from the initial 
condition To,« = f\P n (x)P n (-x)dx = 2(-l)"/(2re + l),n e Z>o, and build Eq. ll9| m ^il inductively from the recursion rela- 
tion given by Eq.fTSl 
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The interpolation formula (known as Dougall's expansion in [Ref.lia p. 167]) 



Mx) = 2 ( 
Jo 



^0 H (cos^-x)cos^2liiM d p 



o \/2(cos/3-x) "■ jo £ = o 



11=2(1 

71 JO f =f 



^ oo (2£+l)/3 (2v+l)^d J 6 

P((X)COS COS 



e=o 



sin(^-v)7r sin(^ + v+ 1)k 1 
+ ■ 



(£-v)n 



(£ + v + l)n 



allows us to sum Eq. |19| m yil as 

/l oo 
-i ^=0 



sin(^ - yL)n ^ sin(£ + /i + l)n 



for [i e C \ Z,n e Z>o. To identify the last term in the equation above with the expression in Eq. ll9| ; , we note that the 
sum over £ in fact truncates after finite terms, which reveals the expression 



sin(pjr) 



(-If- 



r(i^)] z [ r (H±2)] r(^H)r(2^) ^ 



sin(£ - fi)ji sin(£ + \i + l)n 



(£ + H+l)n 



as a rational function of fi, for whatever integer n. As [i approaches any integer m, it is clear that such a rational 
function tends to zero with 0(/i - m) convergence rate, so it must be identically vanishing. 

By Eq. [11] along with the Mehler-Dirichlet theory, we have the following computation when v is not an integer: 



r2n d( k 

P v (x)P v (-x)= I P v ((l-x 2 )cos(p-x 2 )-^- 

JO 27T 

r 2n ' 



I f 

Jo Jo 

Jr*Z7l p7l OO 
/ V P ? ((l-x 2 )cos0-a: 2 )cos 
o Jo f =n 



* flffCcosyS-U-je*) cos (p + x*)cas^Y^ d/3 
v / 2(cos/3-(l-x 2 )cos0 + x 2 ) n 



d0 
~2n 



(2£ + l)B (2v+l)/3d/3 
cos 



d</> 
2n 



= £ Pe(x)PA-x) 

e=o 



sin(^-v);r sin(^ + v+l)7r 
+ 



(£-v)7t (£ + v+l)n 

where we have integrated over ft and <p separately in the last step. Thus, we can perform a decomposition 

sinCvjrXT^v + _ v -i) 



e=o 



sin(^-v)jr sin(^ + v+l)7T 
+ 



(£-v)n 



(£ + v+l)n 



where 



-p£*)r( 



2 J 1 2sin(/i7r) 



W 5 or j2£^- fi j r j2£^±fi)A-/ A* 2 



1-fi 



fi + 2 



4F3 



1 IzE Ell — \ 
-l, 2 > 2 ' 

2—\i [j+3 * 

V —n~, — o - , -L ~~ A 



follows directly from the definition of 4F3. This verifies Eg. 1191 

We shall prove Eg. ll9|9m7^1 bv a contour integral method [Ref.ll9l p. 105], before going back to Eq. ll9n In the next 
few pages, we shall repeatedly use a convenient notation: for x < x',y > 0, the symbol C y x , represents a counterclockwise 
rectangular contour with vertices x — iy, x' — iy, x' + iy and x + iy. 

Akin to the foregoing interpolation arguments, we have 

sin(v7r) 52 n 
T2m,v- }_ (-1) T% m n 

71 



n—m 

sin(v7r) 



v-n v+n+1 

2n+l-2m 



E 



)T(n+l + m) I 1 



[r [i=M)] 2 [T(m+l)f n=m T(n + 1 - m)T ( gn±|±gzz» ) [v-n v + n+1 



v e C. 



Without loss of generality, we may first assume that v e (m, +00) \ Z. For any positive integer N > m, we may establish 
the identity 



r( 



2z + l-2m 



fc N , . T(Z + l-mW{ 



) T(z + 1 + m) n COt(7Tz) dz 71 COt(V7T)r ( 



l-2m+2v 



2z+3+2m 



J v-z 2ni 



r(v+i-m)r( 



-r, — )r(v+i+m) ^ r( 
- + E 



2n+l-2m 



)r(re+l + m) 1 



2m+2v+3l 



n—m 



r(n + i-m)r( 2 " +3 2 +2m ) v-i» 
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by picking up the residues at z = v and z £ Z n [to, AH- Consequently, we have 



>Jc j 



r( 



2z + l-2m 



lim _ 

iv-oo/c^ , . T(z+l-m)T{ 

■ l-2m+2v\ 



)r(z + l + 7n)jrcot(7rz) dz 



2z+3+2m 



j v-z 2ni 



- lim 

iV-oo 



rm-± + 
Jm—i—i. 



iN 



r( 



2z + l-2m 



) T(z + 1 + to) ncot(nz) dz 



N T(z + l-m)r{ 



2z+3+2m 



j v — z 2ni 



7rcot(v^)r ( 1 ~ 2 ^ +2v ) r(v + 1 + m) °° r( 2 " +1 2 - 2m )r(» + i + TO) 1 
r(v+i-TO)r( 2 ^±2) + ^ 



- 2 , n - m nn+l-m)T{ 2n+ l +2m ) v-n 



(20a) 



as the contributions from the other three sides of the rectangular contour vanish asymptotically when the positive 
integer TV goes to infinity, due to the following bound estimates on the integrand: 



r( 



2z + l-2m 



)T(z + l + m) 



r(z + l-m)r( 



2z+3+2m 



= — + 

\z\ 



1 



as z — ► oo; 



(2n + l + 2iy)ji 
cot — 



cosh(2«7r) + cos(27rx) 

= tanh(7ry) < 1, |cot(7rx+ inn)\ = \ <coth7r, Vx,y £ R,n £ Z \ {0}. 

cosh(2/i7r) - cos(27rx) 



2 

Likewise, by a limit procedure, we have 

r( 2 ' +1 - 2 " )r(z+l + m)7rcot(7rz) dz r-m-^iN T ( 2z+1 " 2m ) T(z + 1 + to) yreotfaz) dz 

hm 



,. / 1 r 2 Jnz+l + TO)7rcot(7rz) dz , /•-"»- f+l 
hm <J> r — t—z = hm / 

^ocJcN r(z+l-TO)r( 22+3 +2 ' n ) V+1 + Z27TZ N^ooJ- m -i-it 

r( 2ra+1 2 - 2m )r(» + i + TO) i 



./r:«v , . r(z+l-TO)r( 2z+3 9 +2,n ) v+1 + z2t« N^J- m -i-iN T(z + l-m)T( 2z+s 9 +2m ) v+l + z2ni 

(2n + \-2m\ 



n=mT(n+l-m)T{ 2n+ l +2m ) v+l + n 
where the pole z = -v - 1 is no longer inside the contour C_ m _y 2 w+i/2 - Adding up Eqs. l20al and l20bl we may deduce 
7rcot(v7r)r ( 1 " 2 ^ +2v ) T(v + 1 + to) °° T{ 2n+ \- 2m )T(n + l + m) I 1 1 



(20b) 



■ + 



E 



r(v+l-TO)r( 2m+2v+3 ) n =mT(n + l-m)T{ 2n+ % +2m ) \v-n v+l + nj 

"°° r(i + 2TO + iy)T(iy) tanh(Try) . f°° T{\ + iy)T{iy -2m) tanh(jry) 



r°° H§ + 2m + iyMUy) tanh(Try) f°° l( 

i i — =-= dy+f / — 

J-oo Hi + i v)r(l + 2to + i v) 2(to - v) - 1 + 2 £y J-oo r(± 



-oo T(j + iy)T(l + 2m + iy) 2(to - v) - 1 + 2iy J-^ r(§ -2m + iy)T(l + iy) 2(to - v) - 1 - 2iy 
f°° T(i + 2TO + iy)r(iy) tanh(Tiy) . f°° T(\- iy)T(-iy -2m) tanh(jry) 



— dy 



r°° Hf + 2TO + iy)l(iy) tanhUy) f 00 Kf 

s / — ; dy-i / — q- 

J-oo T(i + zy)r(l + 2to + iy) 2(to - v) - 1 + 2iy J-oo f(i 



T(| + ty)r(l + 2to + iy) 2(to - v) - 1 + 2iy J-oo I{\-2m -iy)T(l - iy) 2(to - v) - 1 + 2iy 
where we have used a consequence of Euler's reflection formula T(z)T(l -z)- 7r/sin(jrz), namely, 

r(i+2TO + iy)r(iy) T(±-iy)T(-iy-2m) 



dy = 0, 



T(i + iy)T(l + 2m + iy) T(\-2m- iy)T(l -iy)' 
in the last step0 Thus, the relation 



TO £ Z, 



7rcos(v;r)r( 2v+1 - 2,K ) T(v + 1 + to) 

T2mv- 5 , VE(TO,+Oo)\Z 

[r (±=fs)] 2 [r(TO+ 1)] 2 r(v + 1 - m)r (^f 2 ^) 

is verified, and one may confirm Eg. H9|9„T!/il bv analytic continuation in v. Once Eg. I19)9m~vi1 is given, and T2 m +i,v = 
0,to £ Z>o is obvious from symmetry considerations, we can interpolate T^ v as follows: 



sin(£ - n)n sin(£ + + l)n 



{£-[i)n (\e + n+l)n 



sin(/i7r) 



' v V - 2m 2m + /u + 1 



r( 2y+1 - 2m )r(v+i + TO) / i 

= sin(/iK)cos(v7r) 2_ 



to [r (l^a)] 2 [Hto+ l)] 2 r(v + 1 - m)T ( 2v+3 + 2m ) \»-2m 2m + fi+l 



^As the cancellation result here is sensitive to the integer constraint on m, the contour integral computation for T2 m ,v, m £ Z>o does not readily 
adapt to the reduction of T^ v with arbitrary /j e C. 
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which in turn, directly accounts for the hypergeometric summation in E q . 1 19 f I concerning 5P4. Exploiting again Euler's 
reflection formula T(z)r(l - z) = Ji/sm(nz), one can rewrite the expression for T2 m ,v as (taking appropriate limits when 
the resulting fractions assume indeterminate forms) 



T2m,v : 



sin(vTT) [T{m+\)fT(m-v)T(m + v + l) 
n [T(m+l)fT(m-v+hnm + v+h , 



m e Z>o, 



which specializes to Eq. H9| 2;u ~ \ i^l 

To prove Eq. I19|„ „J we show that the function 



X(z) := 



sin (nz) 



l + 2cos(«) 7rr(^i)r(^t2) 



[r(¥)l [r(f irrm 

-9/4 



:1(-Z-1) 



(21) 



3 rwl-zll 2 rWz+2tl 3 W 32+3 1 

is analytic over the whole complex z-plane, with asymptotic behavior X(z) = 0(iV~ 9/4 ),z £ C^_ 4Ar)/4 ( 1+4N y 4 as the positive 



integer N goes to infinity. 

First, to verify that T(z) is an entire function, we only need to check that the expression inside the braces of Eq.l2Tl 
has vanishing derivative whenever z is an integer, so that the numerator of Eq. [21] encounters (at least) a second-order 
zero at every integer z = n £ Z. At the positive odd integers z = 2n + l,n £ Z>o, one may directly compute 



dT z 



dz 



= 2 / 

z=2n+l J-l 



1 dP z (x) 



dz 



P2n+l(x)P2n+l(- 



z=2n+\ 



-x)dx + j [P2n- 



2 dP z (x) 

i(-x)i 



dz 



dT Zt 2n + l 



dx = 

z=2n+l OZ 



= 



z=2n+\ 



from Eq.119 ^, „H which coincides with the behavior 

l + 2cos(;rz) nT[£&)T{&&) 

3 [r(¥)] [r(f )] I'W) 

attributed to the factor [r((l - z)/2)] 2 = 0((z - (2n + 1))~ 2 ). Meanwhile, at the non-negative even integers z = 2n,n £ Z>o, 



dT, 



dz 



/1 dp (x) r 1 
— ^ P 2 „(x)P 2n (-x)dx + [P 2 „(-x)J 
. „. 1 z=2n J-l 



2 5P 2 (X) 



3z 



dx - 3 — 

z=2n OZ 



z=2n 



OZ 



log- 



*= 2 » _ [r(^)] 2 [r(^)] 2 r(^)r(^) 



is also a result of Eq. [T9"L „J As we have 



log- 



r(*^=a)r(3J^±£) 



5 r(2±i)r(^f2) V7r 

'"f, , log 5 5 = -27rtan — , 

™ [r(¥)l [r(^)] 2 r(^)r(^-) <5z 2=v ir^firffjfrff ] 2 



which vanishes when v = 2n is a non-negative even number, it is clear that the numerator of Eq. [21] has 0((z - 2n) 2 ) 
behavior for n £ Z>o. 

Then, we investigate the asymptotic behavior of X(z) as z — oo. For any point z sitting on the right half of the contour 
C(i-4A0/4 (1+4A0/4 ( wnere N is a large positive integer) satisfying Rez > -\, we have a uniform asymptotic formula 



Pz,z — 



1 + 



°£)]f 



Jo z + - fi [1+ cos(7rz)] + Yq z + - sin(7rz) YJq \\z + - 



Jq z + - cos(7rz) + Yq z + 



21 ) 



sin(jrz) 



3 



ae. 



(22) 



Here, we have used the following transformations for —n < argv < n: 

-ji/2 



v,v — / 

Jo 

-!" 

Jo 



(cos0)fP v (-cos0)sin0d0 



r>n/z 

[P, 

JO 



(cos0) + P v (-cos0)]P v (cos0)P v (-cos0)sin0d0 



P v (cos 0)[1 + cos(v7r)] Q v (cos0) sin(vsr) }- P v (cos 0) 

jr 



P v (cos0)cos(v7r) Qv(cos0)sin(v7r) 

71 



sin0d0, 
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along with the asymptotic formulae (see [Ref.Efll Chap. 12, Eqs. 12.18 and 12.25] or [Ref.EJ Eqs. 43 and 46]): 

1 



P V (COS0): 



Jo v+-fl+0 



2v+l 



n I 6 
, Q v (cos0) = --i 



Yo\\v + -\e\ + 



2v+l 



2 V sinf 

where the error bounds are uniform for the angular variables 6 e (0,7r/2], so long as |v| — ► oo,-n < argv < n §,H. For 
the case with degree v = (1 + 4A0/4 where AT is a large positive integer, we split the integral on the right-hand side of 
Eq.l22"linto two parts: 



pn/2 pl/V2v+l pn/2 

(•••)d0 = (•••)d0 + 

Jo Jo Ji/v^+T 



(•••)dfl. 



For the first portion, we may deduce 



I 



W2v+1 



Jo v+- \6 [1 + cos(vtt)] + F v+- \sm(vji)\J v+- 



Jq v + - cos(v7r) + Yq v + - sin(vjr) 



3 , „ 8[1 + cos(v7r)][l + 2 cos(vjt)] 
d0 = - 



sin0 



3V3n(l + 2v) 2 



1 + 



from three integral formulae valid for a > [Ref.El items 2.12.42.25, 2.13.22.11, 2.13.25.1]@ 



Jo 



[Jo(ax)] xdx = 



V3na 2 



f 

Jo 



> POO 

[Jo(ax)] 2 Yo(ax)xdx — 0, / Jo(ax)[Yo(ax)] 2 xd;e : 

Jo 



3v/3jra 2 ' 



as well as the familiar asymptotic behavior of Bessel functions for large arguments x » 1: 



which allows us to verify 



/ 

Ji 



Jo(x)~ \J^cos\x- , Yo(x) ~ y^sin|x- . 



Jo v+- 1 [1 + cos(v7r)] + Yq v + - x sin(v7r) > Jo v + - \x x 



x f Jo f[ v + t:\ x \ cos(v7r) + Yq ( (v + -)x) sin(v7r) 



2J J 



xdx = O 



(2v+l) 9/4 



via integration by parts0 Meanwhile, using the asymptotic forms of Bessel functions and integration by parts, one can 
also verify tha10 



pji/2 

JlA/2v+l 



J v + - [1 + cos(vtt)] + Yq V + - sin(VTT) J v + 



x Jo v + - cos(v7r) + Yq v + - sin(vjr) 



3 



■d0 



smt 



3/2 r n-a/V2^+l) 



tt(2v+1) 



■> /■71-U/v: 



COS V + 



cos v + - mr-0) 



tt\ d0 







4j^shi0 l(2v+l) 9/4 



The first integral here belongs to the generalized Weber-Schafheitlin type [Ref. 

El 

§13.46]. The proof of the second integral can be found 
in [23]. The third integral can be derived from the first one, in conjunction with Nicholson's formula [Ref. 15, §13.73] [Jq(x)] 2 + [YoM] 2 = 
-ny Jq° Ko(2xsinht)dt and a special case of the modified Weber-Schafheitlin integral [Ref. Ha. §13.45] Jq(x)K o(2x sinh i)x dx = (1 + 4sinh 2 1) . 



3 To wit, we need to estimate 

f°° sin(a(2v + !)* + &) 



f°° sin(a(2v + l)x + 6> d _ sin 6 f 
Jl/V2v+T \fx aV2v + lJ\ 

cos(av'2v + 1 + 6) sinfe f°° sin(aj>) 



sinfe r°° 15 



- , sinCayJdj' _ 

2v+l sfy oy aV2 



cos b f°° 15 



sb r 
vTliv 



2v+i sfy dy 



cos(av)dy 



b r° 
7TTU: 



Ay- 



cos6 



cos(ay) 



&y = 



1 



(2v + 1) 



3/4 ' 



a(2v + l) 3/4 2a\/2v+ 1 JV2v+T ;y 3/2 2av / 2vTT J\/2v+T y 3/2 

where the coefficients a, 6 e H arise from the product-to-sum formulae of trigonometric functions applied to asymptotic expansions for the products of 
Bessel functions. 

4 One might wish to compare this with the direct application of a hypergeometric asymptotic formula for Legendre functions (see [Ref. 8, p. 295] 
or [Ref. 18, p. 162]): 



P v (cos0) = 



r(v + l) 



r(v+|) v nsmt 



1\ „ n 



where fle[£,i-c],£> 0. 
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Thus, the asymptotic expansion 



8[1 + cos(v7r)][l + 2cos(v7t)] 
3V3n(l + 2v) 2 



1 + 



is established for v = (1 + 4A0/4 > 0. For an arbitrary point z satisfying Rez > -1/2, we may modify the aforementioned 
procedure by deforming the integral path joining the points and n/2 into the sum of two line segments: 



rn/2 r 

(•••)d0 = 
Jo Jo 



v/|2z + l|/(2z + l) (-7T/2 

(•••)d0 + 

J v 



Vl2z + l|/(2z + l) 



(-)de, 



followed by the integrations f£°(---)dx and f^z+H/(2z+i)^" a ^ on S contours that run to infinity in such a manner 
that (2z + l)x — ► +00 along the positive real axis. This results in 



8[1 + cos(;rz)][l + 2cos(7rz)] 
3V3n(l + 2z) 2 



1 + 



Z E C(l-4A0/4,(l+4A0/4'R e2: > 9 ' 



(23a) 



where the error bound is uniform along the contour CjJ_ w)/4 h+anva' on wn i cn D °th I tan(;rz)| and | cot(7rz)| are uniformly 
bounded. After reflection z >-► —z - 1, we obtain 



8[1 - cos(7rz)][l - 2cos(7rz)] 



1 + 



3v^Jr(l + 22) 2 

We may combine Eqs. l23a land l23bl into a single formula 
m l + 2cos0rz) nl[^)T[^) 



zeC 



N Rez<-- 



[r^lflrfflfrlf) 



zeC 



N 

(l-4A0/4,(l+4iV)/4> 



(23b) 



which also implies the bound estimate X(z) = 0(N ) uniformly applicable to all the points z on the rectangular 
contour C^_ 4N y 4 (1+4N y 4 , thanks to the inequalities 



sup 



I cot(7rz)| < max(l,coth7r) = coth7r, 



sup 



zeC 



J (l-4A0/4,(l+4iV)/4 * ^(l-4AT)/4,(l+4iV)/4 

for N e Z>o- By Cauchy's integral formula, we have the identity 

%(z) dz 



I sin(7rz) 



< max 



V2, 



sinh7T 



= V2, 



T(v)= lim <p 

N-coJc 



c. N z — v 2ni 

^U-4iV)/4,(l+4iV)/4 



= 0, VveC, 



which proves Eq. |19|„ vTl 

Both Eqs. |19|9v-i v^l and |19)9,/^9 vTl follow directly from the recursion relation (Eq. H5) . ■ 

Remark As in our proof of Eq. I19|, V ,J we can use Bessel functions to establish the following asymptotic behavior for 
\z\ -^oo,Rez> -1/2: 



8 cos(p7rz + nz) + cos(7rz) + 16 sin(p7rz)sin(7rz) 
PZ ' Z n pv^^2(l + 2z) 2 J? 2 p v / 4 r p f d + 2z) 2 



n - arcsin \ / 1 - — 



8 sin(pjrz-jrz)-sin(7rz) 16 sin(p7rz)sin(7rz) . p 

sinh 



■ pz,z 



" p^P 2 - 4(1 + 2z) 2 p^P 2 - 4(1 + 2z) 2 



0<p<2; 
p>2. 
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along with an asymptotic reflection formula T p( ^_ z _ 1) _ z _ 1 ~ e" t(1 " p)arg2/|arg2| T p2 . 2 .,0 < | argz| < n for p > 0, \z\ — ■ oo. Put 
differently, for large |z| and —n < argz < n, we have an "asymptotic trigonometric modulation" of the standard Clebsch- 
Gordan integral formulae: 



■ pZ,Z 



- pz,z 



r(l^£)]'[r(^)] r(H±<^£)r(<^^) 

cos(p7rz + ;rz) + cos(7rz) 2 sin(p7rz)sin(;rz) 

l + COS(p7Tz) 71 1 + COS(p7Tz) 

wr (i±<^)r(<£±|2±2) 
'(^)] 2 [r(^)] 2 r(^)r(<^) X 



7r - arcsin \ / 1 



0<p<2; 



sm(pjiz-7iz)-sm(jiz) 2 sin(p7rz)sin(7rz) _-, / p ; 

smh \ 1 

4 



l + COs(p7Tz) 



71 1 + COS(p7Tz) 



cot 



(p-2)nz 



p>2. 



For p ^ 1, the right-hand side of the penultimate (resp. last) asymptotic formula diverges at z = -2/(2 + p) (resp. z = 
-Hp), so the asymptotic analysis does not result in exact forms of T pZtZ (p ^ D for finitely sized |z|. □ 



Corollary 2.2 (Some Multiple Elliptic Integrals in Clebsch-Gordan Forms) We have the integral formulae 



— T 



^-T _i/2 = 2 f 1 K(^^t)K(^/^~t)at= 
4 Jo 4 
,2 r i 

— r 01 / 2 = 2 [2E(v / f)-K(v / f)][2E(v / 1^7)-K(v / l :: 7)]di = 
4 ' Jo 

/■l _ [r(i)] 8 

mVT r t)] 2 K(Vi)dt=—^ T , 

Jo 1927r z 
27 f 1 (l-p 2 )p(2 + p) [ / /p 3 (2 + p)^ 



■1/2,-1/2 - 



(24) 
(25) 
(26) 



-77 T -113 



-1/3 : 



-T- 1/4- 1/4 = 



27 r 1 

2 Jo ^3 + 6p(l + p+p 2 ) 

= 27 f 1 (l-p 2 )p(2 + p) 
6 Jo ^/T+2^(l + p+p 2 ) 

r 

Jo ft -l ,/^|3/2 



K 



K 



l + 2p j 



P 3 (2 + p)\ 





2 






K 











p 3 (2 + p)\ 



-1/6,-1/6 - 



(1+ Vu) 3 
Jo ll + v^/ 

£ 



K 



J 1 + Vu 



K 



l + 2p 

2v^ 1 



K 



l + 2p j 
p 3 (2 + p)\ 



l + 2p j 



dp 
dp = 



256?r 2 ' 



(27) 



dw = 4v / 2 



(l + f) 3/2 



At 



K 



2-Ju 



27 
2^ 



1-x 



1 + 

2 / 
K 



K 



^ V 1+ 



lu = 4 f 
Jo 



Mi - ^)[K(vT^7)] 2 K(y^) [r(i )r(|)] 2 
o U + *) 3/2 24 ' 



l + x 



1-x 



(3 + 



_ 27 r 

^4 TJo 



27 f 1 t(l-t)[K(VT^t)] 2 K(y/i) 



(l-t + t 2 ) V4 



dt-- 



[r(i)] 4 



8^2^ 



-77 T 112 -31 A 
o 

8 



^ r 1 2E(y^)-K( v^) 
Jo 1 + \/w 



-K 



2y^I W / l-y/u 
1 + x/wj \V l + \/w 



du = \/27r, 



2^ /' 1 8(1-2 M )E(v / w')-(5-8m)K(v / ^) 



■ 3/2,-1/4 : 



Jo 



1 + \fu 



K 



2v^ \ 
1 + Vu 



K 



J 1 + Vu, 



du- 



\/2~7l 



(28) 
(29) 
(30) 
(31) 



Proof Special cases of Eas. [19g^7^1ll9^^l and ll9|9v..9 J lead to Eqs.[2l[25l|26l[27j|2l|29l[30]and|3Tl ■ 

3 Spherical Rotations and Multiple Elliptic Integrals 

In this section, we shall focus on the transformations of multiple elliptic integrals using variable substitutions 
motivated by rotations on a unit sphere. In effect, our attention will be momentarily restricted to two special Legendre 
functions P-1/2 and P-1/4. 
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The transformation methods in this section not only provide evaluations of certain Cauchy principal values: 



rilMwjgW _ 1<x<1; 9 fJ 
J-i jrCx-f)v/l+? J-i [ 



1 + t 



I 1 



K 



2d£ 

- 1 < 3c < 1 etc. 



7r(x - £ ) 

which will be used in {j4j but also pave way for further applications in our subsequent works on elliptic integrals. 



3.1 Beltrami Rotations 

The technique of spherical rotation, which traces back to Beltrami's work on the Abel integral equations [Ref. 
p. 328], has found applications in the geometrically-motivated evaluation of certain definite integrals related to Bessel 
functions (see [Ref. 15, §3.33, §12.12 and §12.14]). The lemma below is a simple realization of the Beltrami rotation in 
the case of the complete elliptic integrals K(k) and E(&). 

Lemma 3.1 (Beltrami Transformations) (a) We have the following integral identities 



n Jo l-k A K A 



K(0- 
K(r): 



2 r 1 Ji-( 2 K(VT^)dK 
- r 1, — 2^—> o<f<i, dH) 

ttJo 1-<^(1-jc 2 ) 



" ,/ \ 2 , 2 ■ 0<r<l, 1HL) 

7tJo (1 + rr -4rK z 

m)=- ' >a , a > o<ij<i. 4HL) 

7rJo (l-7?r + 4rjK 2 L 



(b) for < & < 1, we Ziaue 



K(/e)-E(/e) 2 f 1 Etv'T^dje 



_ 2 r+ECV: 
Jo 1- 



(33) 



k 2 71 Jo l-k 2 K 2 

4(1-^ 2 ) r 1 E(v / 1^2)d ) c 
E(/j) 4(1- A 2 ) f 1 E(vT :: K2)d)c 



-g) f 1 E(Vl-^)d* 

?r Jo [1-& 2 (1-jc 2 )] 2 ' ^ 



Proof (a) We start by writing the complete elliptic integral of the first kind as an integral on the unit sphere S 2 : 

d8 1 f n . „ ,„ r 2n , l r 1 do- 



K(A) 



If 71 d0 1 f" , 1 f 

= - ■ — / sin0d0 d0 = 

2 Jo v/l-jfe 2 sin 2 4:71 J° -> sin 0(1 -ft sin cos 0) Js 2 \/l-Z 2 (l-&X) 4?r 



As we rotate about the Y-axis by a right angle, effectively swapping the roles of the X- and Z-axes, we may rewrite the 
integral above as 



f 1 do- f 1 da 1 f 1 r 2n 

K(A)= / ■ = ■ = — sin0d0 d0— 

Js 2 Vl-X 2 (l-kZ) 4^r Js 2 v / l-X 2 (l-ft 2 Z 2 ) 4?r 4;r J Jo 



sin 2 B cos 2 (p{l-k 2 cos 2 0) 

where we have combined the values of the integrand at points +Z and restored the spherical coordinates. Integrating 
over the azimuthal angle 0, we obtain 

- , w K(sin0)sin0d0 

K(k) ■ 



1 r*K(si 
ti Jo 1- 



- k 2 cos 2 
which is equivalent to Eq.[32"l 

Now that we have proved the relation 

m) -A( l ^^l^, <k<i, 

71 Jo l~k z K z 

which analytically continues to k e C\[l,+oo), we may derive Eqs. 13211 132H and 1321 1 respectively, from the imagi- 



nary modulus transformation K(£) = K(if/\/ 1 — % 2 )/V 1 — f 2 , Landen's transformation K(r) = K.(2y/r/(l + r))/(l+ r), and a 
combination thereof K(?y) = K(2iy/rj/(l - J7))/(l - ??)• 
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(b) We adapt our derivation of part [(a)] to accommodate for the complete elliptic integrals of the second kind E(£),0 < 
k<l: 



K(k)-E(k) _ 1 f* sin 2 0d0 _ 1 f" s - n9de f 2 " d sin9 _f v'l-Z 2 da 

k 2 2 Jo \fi-k 2 s\n 2 Q 471 Jo Jo l-£sin0cos</> Js 2 1-kX 4n 

Vl-X 2 da _ f Vl-X 2 da _ 1 f n _._ aJO f 2n j t \/l - sin 2 0cos 2 </> _ 1 f M E(sin0)sin0d0 

£ 2 COS 2 



r Vl-A z d(T_r Vl-A z d(T_ 1 /" in0d0 Id V1_sin 0COS ^ - 1 f E(si 
"Js 2 1-&Z 4n~Js 2 l-k 2 Z 2 ^n~4nJo Sm Jo l-/z 2 cos 2 ~ Wo ~ T 



which verifies Eq.[33] As we multiply both sides of Eq.[33]by k 2 , and differentiate in k, we obtain the formula stated in 
Eq.m 



AE(A) d rT „„ 2d f 1 E(vT^K 2 )^ 2 dK 4* r 1 B(Vl^it 3 ')dK „ . , 

-ITO)-Kft)] = --J o =-J o a _ k2x2)2 , 0<k<l. 



1-k 2 dk 
Now, by analytic continuation, we have 



4(l-& 2 ) f 1 E(vT^ 2 )d)c 

E(/s) = 



-k<) f'-EWl-K^dK 

~ Jo (1-^K 2 ) 2 ' ^^^^ ^ 



In particular, by the imaginary modulus transformation for complete elliptic integrals of the second kind [Ref. |2£ 
item 160.02], we obtain 



4(1 + ^ 2 ) /■ 1 E(v / l-K 2 )dK 



l^ITf 2 / « Jo (1 + < 2 k 2 ) 2 



which is equivalent to Eq.|34 



Remark Admittedly, the spherical rotation is not the only road to any or all of the identities in Eqs. [32j 13211 132H and 
132^1 For example, one may still verify Eq.[32]via the method of moments and hypergeometric summations [6]. The same 
can be said for several formulae that we will encounter in Proposition 13.31 However, it is our hope that the spherical 
rotations provide clearer geometric motivations than a purely combinatorial approach. □ 

Beltrami transformations allow us to convert one multiple elliptic integral into another. For instance, the following 
chain of identities 



fmwf''f a d..f' K( "L_ ^.'fWia 

Jo Jo 1 + Jc Jo 1+ Jl-E 2 Jl-E 2 nJo \ ) 



yi^Fx/w 2 " 

result from two applications of Eq. [3211 and the transformations 



JcarccosK , 

-log(2y) 

Vl-K 2 



dx 



f 1 ™,^, 2 f 1 E(Vl-K 2 ) r , 2n f 1 (2 + ?c)E(v / l-Jc 2 ) n 

/ E(k)dk = -l — — 5 -[-K + (l + ?c 2 )tanh 1 x]dx= / -—^5 -dx 

Jo JtJo f Jo (1 + kY 



can be justified by Eos.l34land l34ll More examples will be described elsewhere. 

In the next corollary, we mention a few consequences of Beltrami transformations that are pertinent to later devel- 
opments in the current work. 

Corollary 3.2 (Some Applications of Beltrami Transformations) (a) For < u < 1, we have the following duality 
relations for multiple elliptic integrals: 

K(b.Uh ri kK(k)dk 

(35) 



r^ 1 K(k)dk r 1 
Jo VT^Vu^k 2 Jo VT^W^i-k 2 u 
r 1 kK(k)dk _ r 1 K(Vi-K 2 )dK 

— U Jo Vl — K 2 Vl — K 2 U 



(36) 



(b) We have the following Cauchy principal values for - 1 < x < 1: 



9 fi Kiva^m) Kiyft+jm 

J-in(x-!;)y/T+{ y/T+x 

, r 1 K(v/(l-0/2)-2E(v/(l-0/2) Jr _ 2E(v / (l + x)/2)-K(v / (l + x)/2) 

/ , di, . (3a) 

J-l 7l(x-0\/l + l VT+x 
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Proof (a) With Eg. 13211 we may interchange the order of integrations to compute 
K(k)dk f^\2 f 1 x/l-$ 2 K(VT^K2)dK 



r vu K(k)dk _ n u 2 r 1 y/i-j 

Jo Vl-k 2 Vu-k 2 Jo nJo 1- 



<u _ r 1 K(v / r r ??)djc 

s/X-ps/u-i 2 Jo sJl-u{l-K 2 ) 



where the last expression is equivalent to the right-hand side of Eq. [35] by the correspondence k Vl-k 2 . Similarly, 
we may prove Eq.[36]using Eq.1321 



f 



kK(k)dk 



_ r v ^[2 r 1 K(v / l-)c 2 )dK 
w Jo ttJo 1-& 2 )c 2 



==( 

■u Jo 



1 K(v / 1-K 2 )d)f 



v / 1-7c 2 v / 1-k: 2 u 



(b) For z > 1, one can readily compute (cf. [Ref.ll9l p. 233]) 

dt r 1 '^ dt 



Jo 



-f 

Jo 



JllsJ 



dt 



Vl-t 2 Vl-zt 2 Jo Vl-t 2 Vl-zt 2 JVs/i Vl-t 2 Vzt 2 -l 
1 ds i fv^ ds K(VT/z)-iKW(z-D/z) 



1 r l ds i n z ds 



This is the inverse modulus transformation for complete elliptic integrals of the first kind. 
Now, in the identity 



(39) 



K(k): 



2 r 1 K(VT^)dK uu 
-\ — ; — To~^ — > VfeeC\[l,+oo) 

n Jo l-ft z K z 



we approach the limit of k + i0 + E [l,+oo), read off the real part, and employ the inverse modulus transformation 
(Eq.MB, to obtain 

k r 1 K(Vl-x 2 )dK 

K(k)=-@> — —„ i- — , 0</e<l, 

n J-\ 



k 2 -K 2 



which can be reformulated into Eq. [37] after the variable substitutions k = \/{l + x)l2, k = \J{1 + Z,)I2. 
From Eq.[37l we may directly compute (cf. [Ref.S Eq. 11.215]) 



1 



1 K(v/(l-0/2)VT+Td£ _ f 1 K(v/(l-0/2)[l + (f-x) + x]df 



= (1 + 



f 1 K(V(l-Q/2)d< H 1 K(V(l-Q/2 )d^ 

X i-i rc(x-0 v / i+? fi-i \/T+? 
rr+T| l r 1 K(v/(i-0/2)df 



1 < X < 1. 



Here, the last integral can be simplified with the knowledge of / K(VT^K 2 )dK = n 2 /4 \Ret. Ila. item 6.141.2], which 
brings us 



K(V(l-Q/2)^I+?d^ 
?r(x - () 



Vl + xK 



1 + x 



v 2 ; 



— , -i<x<i. 



Applying integration by parts to the Tricomi transform (cf. [Ref. 26, Eqs. 11.218 and 11.219]), one can use the equation 
above to deduce that 



•f- 

J-l 7t( 



1 d 



In other words, we have 



K 



f 1-f 



VI 2 I 



vl + xK 



1 + x 



v^(x-i)' 



-1 < X < 1. 



^ K(v/(l-0/2)-E(v/(l-0/2) Jr 1 

■d< = 



J-l Jp 



(x-oa-ovi+z 1 ~ x 



E(V(l + x)/2) 1 
Vl + x \/2 



-1 < x < 1. 
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Carrying this further, we obtain 

9 L 



1 K( v/(l-Q/2)-E(^(l-0/2) d ^_ E(V(l + x)/2) 1 | 1 r 1 K(y/(l-Q/2)-E(^(l-0/2) 



(i-0\/T+? 



EWd + x)/2) 1 1 f 1 d 



xV2)_J_ 1 r _d_ 
x V2 TiJ-idt 



K 



Jr E(v / d + x)/2) 

= ; , -1<X<1. 



We may combine the relation above with Eq.[37]into a symmetric form 
as stated in Eq.|38l 



1 K( y/(l - fl/2) - 2E( y/q - d ^ _ 2E(y / (l + x)/2)-K(y / (l + x)/2) 

1 7T(X-0\A+T VT+X 



-1 < X < 1, 



3.2 Ramanujan Rotations 

The next proposition is based on reverse engineering of some formulae in Ramanujan's notebooks that were stated 
without proofs. 

Proposition 3.3 (Ramanujan Transformations) (a) For s,t e [0,1), we may represent K(\/s)K(\/i) and [K(\/t)] 2 as 
integrals on the unit sphere S 2 : 



K(v^)K(v^) = i f -= 
[K(n/7)] 2 = f 



da 



v/(l-sX 2 )(l-iY 2 )-Z 2 ' 
2(2-0K(v / X 2 + Y 2 )do- 
ga (2 - £) 2 - t 2 X 2 4^r 
(l + ^KXv^ + Y 2 ) do- 
's* (l + f) 2 -4fX 2 In' 



(40) 
(41) 
(EL) 



where X - sin0cos</>,Y = sin$sin0,Z = cos0 are the Cartesian coordinates on the unit sphere S 2 . Consequently, the 
following identities hold for 0<u<l: 



r^ 1 K(k)dk _ r 1 
Jo VT^k 2 Vu~^k 2 ~ Jo 

rVT^ K(k)dk _ r 1 
Jo \/l-k 2 Vl-U-k 2 Jo VT^WyJ l-k 2 (l -u) 1 + \fu 



kK(k)dk 



VT^W\/i-k 2 u i + Vu 

kK(k)dk 



(b) We have the integral identity 



( / 2v^ 



k V i + Vu, 

K 



J l + s/u, 



(42) 
(43) 



J s K(Vx 2 + Y 2 )f(\X\) < ^-=- I f(k)K 



An 



2 r 1 
= ~ / ' 

71 Jo 



f 1 






1 + k 


K 






m 


[u] 













\dk 



(44) 



so long as the function f(\X\),0 < \X\ < 1 assumes non-negative values and the surface integral is convergent. This allows 
us to convert Eq. ^]]into the following identities for < t < 1: 



[K(v^)] 2 = - [ 1 ^v^v 71 ^^ 

71 Jo \-\lt 



lit 

8 f 1 K(^)K(^T^)d^ 
(l + v^-Ml-v^) 2 ' 



if 

71 JO 



dm 
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and to deduce the following integral formulae^ 



r 1 K(v / T^2) ] [K(VT=t)-iK(Vt)] 2 n 1 

/ (\k = , 0<t<-; (45) 

Jo 



f 

Jo 



V(2t-1) 2 -(1-K 2 ) 2 ' 2' 

K(vT^2) [K(VT~t) + iK(Vt)] 2 

dK= , i<0. 1 1451 ) 



V(2i-D 2 -(1-K 2 ) 
Moreover, we have the identity below for < A < 1/2: 

K(k)dk 



Jl-2A A/l-fe2./fc2_n -2A) 2 Jd- 



-2A vT :: P'v / ^ 2 -(l-2A) 2 J(l-A)/(1+A) vT^^Cl + A) 2 £ 2 - (1 - A) 2 ' 

or equivalently, 

■■>'» K(sin0)d0 _^./T^./^-r,_ r /2 Jfl r /2 <*0 



p K(sm0)d0 = K(v ^ )K(vCT) = p de p 
"0 \/l-(l-2A) 2 cos 2 Jo Jo y^jT 



(46) 



61) 



. + A) 2 - (1 - A) 2 sin 2 - 4A sin 2 < 
Proof (a) Writing da - sin0d0d0 for the surface element on the unit sphere, we have 

-2m a ,h 1 /• (Jg- 

■•(f) o Js2 \/l - sZ Vl - Z 2 - *X 2 

where the Cartesian coordinates X = sin0cos0,Z = cos (9 were used in the last step. As in Lemma I3TT1 we interchange 
the roles of the X- and Z-axes, to deduce 

k(v^4 ( - 7 A ; - r'\e r'\, 

8 Js 2 Vl-sX 2 Vl-X 2 -tZ 2 Jo Jo 



K (v ^) = i r ~jl= r - \ r 

8 Jo vl-scos 2 Jo Vl-fcos 2 ^ 8Js '' 



s 2 v / l-sX 2 v / l-Z 2 -iZ 2 Jo Jo - s sin 2 cos 2 0^1 - sin 2 cos 2 - 1 cos 2 

71/2 sin0d0 r Jl/2 d0 



f W5! sin0d0 
Jo Vl-tcos 2 6 Jo /7 



- s sin 2 cos 2 J 1 - = s l n e ia cos 2 

For any two numbers a, 6 £ (0, 1), we can verify the identity 

-n/2 d rjt/2 d y 



rim A( p ^ rnrz 

Jo J l - acos 2 (bJl-b cos 2 Jo 



yl — a cos 2 0\/l - 6 cos 2 Jo \/l _a_ ( 0_a ) cos2l / / 
with a simple substitution = arctan(\/l-atani//). Now, setting a - ssin 2 0,o = we are led to the formula 

-u/2 rn/2 d0 



(47) 



sin0d0 / —= 
o Jo v(T 



s)(l - sin z cos" 1 0) + (s - ^cos" 1 + s£ cos" 1 sin z sin z 



= if - dg = *f - da 63) 

8 Js 2 y/(i - s)(l - Z 2 ) + (s - i)Z 2 + stY 2 Z 2 8 Js 2 v/(l-sY 2 )(l- £Z 2 )-X 2 

after a literal transformation from the spherical coordinates to Cartesian coordinates, and an algebraic simplification 
according to the spherical constraint X 2 + Y 2 + Z 2 = 1. Clearly, the last term in Eq. l40ll is equivalent to the right-hand 
side of Eq.|40l as the integral in question remains invariant under a cyclic shift of variables (X,Y,Z) >-» (Z,X,Y). 

When s = t, we may swap the X- and Z-axes in the penultimate term of Eq. 1401 1 to derive 

n/2 pn/2 



[K(Vt)?=u da ==r do( n d0— = 

8Js* ^(l- t )(l-Z 2 ) + t 2 Y 2 X 2 Jo Jo v^T7 



+ 1 2 sin 2 sin 2 cos 2 



= -/ d0 / d0 — = -/ 

2 JO Jo / n , , t*-2a „ t „2. k 8j S 2 



l-i+^sin 2 0sin 2 OJS J(l-Z 2 )(l-f+fy 2 ) 



5 Here, in Eg. 1451 we have chosen the univalent branches of the square roots such that Im s/ (2t - l) 2 -(1- jc 2 ) > 0, and K(\/l - f) > 0,K(v7) a 0; in 
Eq. l45ll it is understood that \/(2t - l) 2 - (1 - jc 2 ) > on the left-hand side and 

K(v7)= - >0, ImK(vT :: 7) = Im < 0, Vi<0. 

^0 v/l + U|cos 2 e Jo v/l-d-i)cos 2 e 
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Switching the Y- and Z-axes, we then arrive at 

da 

JS 2 



[K(s/t)] 2 = ~ [ da = fusing d0 f' % d<p— 

8js \/ ( i-Wi-^z 2 ] jo jo xTi 



a-Y 2 )(i-t+%z 2 ) 

r ' 2 K(sin0)sin0d0 f M K(sin0)sin0d0 



-sm 2 6sm 2 (t>\/l-t+ ^cos 2 ' 



r m * K(sin0)sin0d0 _ C 
'° \/l-£+^cos 2 '° 



t+*cos 2 6 Jo V (2 -t) 2 -t 2 sin 2 8 



Judging from the familiar integral formula 

c2n 



d<p 



r M d^ 

Jo 2-t + tsi 



2jt 



it is evident that 



4nJo Jo 2-t 



sin0 cos v^-f^-^sin 2 ^ 
±: < ' ,n 2K(sin0)sin0d0 1 f 2K(VX 2 + Y 2 )da 



t + tsin8cos(f> 

Pairing up the values of the integrand at +X, we obtain 



4n Js 



S 2 2-t + tX 



4% JS 



1 f 2K(VX 2 + Y 2 )da 2-t f 2K(VX 2 + Y 2 )da 

~4~n Js 2 (2-t) 2 -t 2 X 2 



as claimed in Eq. EH Then, Landen's transformation 

K(Vt): 



2Vt 
l + s/t~~\l + s/i 

brings us to Eq. l41H 

On account of Eqs.l4Tland l41hl we may put down the relations 

2 



, 0< t < 1 



(48) 



K 



K 



{ / 2\Hi 



r- r K(VX 2 + Y 2 )da ~ f 1 

d+v^) —i — ^^i-=( 1+v ^W ■ 

Js 2 l-uX z 4ji Jo 



kK(k)dk 



( l-Vu 



t V i + Vu t 



s/l^Vl-k 2 !*' 

l + Vu r K(VX 2 + Y 2 )da f 1 kK(k)dk 



r K( 

Js 2 1- 



(l-u)X 2 4n 



7 

Jo 



VT^W^/i-kHi-u)' 



which confirm, respectively, the right halves of Eqs.[42]and[43] The left halves of Eqs.[42]and[43]follow from the duality 
relation in Eq.[35l 

One may wish to compare Eqs. [41] and I41H (applicable to < u < 1) to the Hobson coupling formula in Eq. 1 13 M 
(confined to < u < 1/2). 



(b) We start with a spherical rotation 

J s K(Vx 2 + Y 2 )f(\X\) ( ^=J s K{Vz 2 + Y 2 )f(\Z\) ( - 

= ^ 1 dZ^ ?r/2 d0K|\/z 2 + (l-Z 2 )sin 2 (/)j f(\Z\)= dk J*'* d(j>Kyk 2 cos 2 <p + sin 2 ^f(k) 



da 
4ji 



=-['«»[[ 

n Jo Jo 



Vl-k 2 -K 2 



dx 



dk=—Re[ f(k) f 
71 Jo Jo 



1 KiVT^K 2 ) 

Vl-k 2 -K 2 



dK 



dk, 



where we have made the substitution k = Vl-k 2 cos(f> in the last line. Then, we may evaluate 

* x K(Vl-K 2 ) 



Re 



/ 

Jo 



JO Vl-k 2 -K 2 

by analytic continuation. From Eg. 1411 we have 

- 1 K(Vl-K 2 ) 



:dK 



-Im 



f 1 KCVT^ k 2 ") 

Jo vprf 



dx 



Jo 



1-Z+^K 2 



'-[/' 

[Jo 



dt 



Vi^WVT^zt 2 



0<z< 1. 
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After setting z — 2/(1 + k) > 1 in the inverse modulus transformation formula (Eq.l39ll. we obtain 



1 r 1 KiVT^x 2 ) l + k 

(l + k) I — Ax 
Jo 



Vk 2 -l + x 2 



K 



l + k 



iK 



1-k 



0<k<l. 



(49) 



Reading off the imaginary parts on both sides of the formula above, we arrive at the following equation for A = (l-k)/2 e 
(0,1/2]: 



kK(k)dk 



1-2A y / l-k 2 \/k 2 -(l-2X) 2 



= K(v^)K(vT^A), 



(50) 



as well as the identity claimed in Eq. l44l 

Applying Eq.[44]to Eq.HH we obtain the result stated in Eq. l41H 



[K(V7)] 2 = ( 
Js 



2(2-t)K(VX 2 + Y 2 )da 
(2-t) 2 -t 2 X 2 



Jo 2-t- 

I 




dk + 



f 

Jo 



2-t-k't 



V^ 2 2K(v^)K(v^7I) 
2-t + {l-2^)t 



d^+ / 



1 2K( V /^ 7 )K( V /1^7) 
1A/2 2-t-(2n'-l)t 




(An equivalent form of Eq. l4TH has appeared as Eq. 26 in [6], which was derived combinatorially using generalized 
hypergeometric series, instead of the geometric interpretation given here.) One may derive Eq. l41£l from Eq. l4iri and 
Landen's transformation 



KiVl^)- 



1-s/i 

i + Vt"{i + Vt 



Writing t = (l-k)/2 e (0,1/2), we may recast Eq.|49"linto Eq.[45] By analytically continuing Eq.|45]to negative valued 
t < with the aid of the integral representation of the elliptic integral K for complex- valued modulus, we arrive at 
Eq. 14511 Here, the change in the sign from -iK(\/f) in Eq.|45]to +iK(\/f) in Eq. l45ll is worthy of special attention. It is 
ready to appreciate that such a sign change ensures compatibility with the natural choices of the univalent branches 
of the square roots occurring in all the places. 

The leftmost equality in Eq.|46]has been already proved in Eq.[50l and its counterpart in Eq. l46ll follows from Eo.l36l 
To demonstrate the rightmost equality in Eo. 1461 we perform the following computations for < A < 1/2: 



f 



K(k)dk 



a-Ava+A) Vi-k 2 x/(i + W 2 k 2 -(i-A) 2 



rn/2 

Re / 

Jo 



K(sin0)d0 



\/4A-(l + A) 2 cos 2 2sfl 



Re 



K 



{ i(l-Vl) 



2VA 



K 



{ 2VI 



by using an analytic continuation of the leftmost term of Eq. l46ll in the last step. According to the imaginary modulus 
transformation and Landen's transformation, we have 



K 



2VI 



2VX ^fi-v 7 ! 



= VXkW1-A); 



whereas the inverse modulus transformation (Eq. 

.((.1 + VX)^ 



1 + y/X + 

and Landen's transformation lead us to 



Re 



K 



2VX 



2_n 

i + vH' 



K 



I 2VX 



l + v r X l 



= 2^AK(v^). 



Thus, the rightmost equality in Eq.[46]is verified. The connection to its counterpart in Eo. l46ll can be proved as follows: 



r 1 ~K.{k)dk _ r nl2 
jvh vT^FvF^ Jo 

-ti/2 pull 



\/l-(l-u)s\n 2 ( 



l-(l-u)sin 2 



dO 



/"7T/Z p 

= d6 
Jo Jo 



dcp 



\/l-(l-M)sin 2 0-usin 2 
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Here, we have used the variable substitution k = Vu/[l-(l-u) sin 8] to complete the first line in the equation above, 
before spelling out K(\/ u/[l - (1 - u)sin 2 8]) as an integral in in the second line. 
The last equality in Eq. |46l| implies the following evaluation for < u < 1: 



Jnn/2 rn/2 
ae 
Jo Vl-u 



d0 



sin 2 0-(l- u)sin 2 



1-Vu 



K 



1 + \/u ) \ V 1 + s/u t 



2\fu~ 



l + s/u 

so an interchange of the angular variables 8 and brings us back to the u <->■ 1 — u symmetry displayed in Eq. [14] 
One may also wish to compare the formula 



Mr) 



l + Vu 



K 



2V~u 



J 1+Vu, 



K 



}l i+\/uj jvh 



K(k)dk 



Q<u<l 



(51) 



to the corresponding result from Hobson coupling (Eq. ll3n . 
Remark (1) Up to a variable substitution, our intermediate result in part (a) 

-" /2 K(sin0)sin0d0 f 1 K(VT :r i^)dK 



_ r 2 r z K(sin0 sm0d0 f K(v 1 — k )cLk 

J ° \/l-i+^COS 2 J ° Jl-t+^-K 2 



<t<l 



(EH) 



is related to formula 13.0.18 in [27], and item 2.16.5.7 of [28]. However, it is highly probable that the ultimate source 
of this formula is Ramanujan. 

Entry 7(x) in Chapter 17 of Ramanujan's second notebook essentially reads 

2 



ji/2 r n/2 



Jo 



dddip 



V(l-u sin 2 8)(1 - u sin 2 8 sin 2 0) 



_ I r" /2 dy/ 
r° v 1 - u sin 4 \fr , 



which was verified by B. C. Berndt with some highly technical transformations of the generalized hypergeometric series 
4^3 [Ref.[Hl pp. 110-111]. In fact, the left-hand side of Ramanujan's formula equals 



1 — u sin 



,/2 K(vw'sin0)d0 _ /" r/2 K(sin0)sin0d0 
~2a in T~Xa 



1 — u sin 



after integration in (p and reference to Eq.[35l while its right-hand side equals 7T 2 [P_i/4(l-2u)] 2 /4, according to the inte- 
gral identity in Eq.[6] In this manner, we see that Ramanujan's formula is actually equivalent to Eq. 14111 a consequence 
of spherical rotations. 

In [Ref. 11, pp. 111-112], Berndt followed up with a combinatorial proof of Entry 7(xi) in Ramanujan's notebook: 



rn/2 rn/2 £sin0d0d0 _ f m C^ink 

J° 'o V(l-k 2 sin 2 0)(1 - k 2 sin 2 8 sin 2 0) ^° ^° 



d8dcf> 



Vl-k 2 sin 2 8 - sin 2 8 sin 2 



K 



1 + / 



K 



1-k 



which turns out to be exactly the real part of our Eq.|45l a geometrically motivated result. 

Berndt has described his proofs of these two entries as "undoubtedly not those found by Ramanujan". Judging from 
the use of spherical coordinates in Ramanujan's presentation, we take leave to think that the technique of spherical 
rotations given in the proof of this proposition might be closer to a rediscovery of the pathway that Ramanujan has 
originally undertaken. 

(2) Combining our analysis in part (b) with Eq. |47l we have effectively shown that spherical geometry entails the 
following integral formula 



r n/2 

K 

JO 



h 2 cos 2 + sin 2 d 



rn/2 

Jo 



K(sin0)d0 

Vl-k 2 cos 2 8 



K 



1 




1 


1 + k 


K 


1-k 










IV - J 
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which is a relation that has nearly 80 years of history. In 2008, Bailey et al. used the method of Bessel moments to 
discover the following integral identity (Eq. 49 in ft) 



I 



71/2 K(sin0)d0 



K 



(sin|)K(cos|), 



V 1 - cos 2 a cos 2 6 * 2> 2> 
and remarked on its equivalence to a formula derived by Glasser in 1976: 

-nl2 



Jo 



K 



sin 2 acos 2 0j dtp — K^sin— jK^cos — j . 



As pointed out by Zucker (30, the left-hand side in the formula of Glasser, being a type of generalized Watson integral, 
can be expressed in terms of 34, an Appell hypergeometric function, which in turn, reduces to the product of two 
complete elliptic integrals of the first kind, according to a result by Bailey in 1933 [29]. □ 

Corollary 3.4 (Some Cauchy Principal Values) For -1 < x < L we have 



'/>( 



1 + f 



K 



2df 



n(x - f ) 



K 



1 + x 



K 



\-x 



(52) 



and 



K 



1 + f 



K 



( l~t\ 



K 



1-x 



+ K 



1-x 



+ K 



E 



1-x 



+ K 



1 + f 



1 + x 



K 



E 



E 



f 1-^ 



df 



n(x - f ) 



1 + x 



(53) 



Proof From Eq.HQand the inverse modulus transformation Re{[K(l/v^)] 2 /i) = \K(y/t)] 2 - [K(VT^t)] 2 for t e (0,1) 
(cf. Eq.[39jl, we may deduce 



[K^)] 2 - [K(v^)] 2 = *V C K( ^ )K( ^ )d ^, 
Jt Jo t - [l 



<t<l, 



which is equivalent to Eq.l52l 

Now, we proceed as in the proof of Corollarv l3.2tb)l and compute 

. , 2 



•/>( 



2(l + f)df 



;r(x - f ) 



= (l + x)< 



:(1 + X)< 



K 



K 



1 + x 



1 + x 



K 



K 



/l 


-x 


2 


/l 


-x 



2 r 1 
n J-i 



K 



1 + f 



K 



J 2 



df 



(54) 



where we have quoted the result T _y 2 = l\P-ii2^)P-ii2^~Od^ = n from Eq. |19|9 m "y>l Likewise, we may deduce from 
the equation above another Cauchy principal value: 



'/>( 



1 + ^ 



K 



2(1 -t )d t=(l-x>)< 



n(x - f ) 



K 



1 + x 



K 



1-x 



7T 2 X 



Like what we did in the proof of Corollary |3.2tb)l we can use the last equation to derive 



d-r)K 



K 



2df _ _d_ 

7r(x-f) dx 



K 



1 + x 



K 



1-x 



7T 2 X 



(55) 



after integration by parts. A little algebra then reveals that Eqs. [54] and [55] together entail Eq. 
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4 Finite Hilbert Transform and Tricomi Pairing 

4.1 Parse val Identity for Tricomi Transforms 

The finite Hilbert transform on the interval (-1,1), also known as the Tricomi transform, is defined via a Cauchy 
principal value (see [Ref. EE Chap. 4] or [Ref.HE Chap. 11]): 



(STf)( 



r f 

x):=&> - 

J-l 71 



1 now 

ix-O' 



a.e. xe(-1,1). 



This induces a continuous linear operator 3~ :L P (-1, 1) — ► L p (-1, 1) for 1 < p < +oo [Ref. l3ll p. 188]. 

We can derive some new integral relations from a Parseval-type identity for the Tricomi transform (see [Ref. [3C 
§4.3, Eq. 2] or [Ref. 26, Eq. 11.237]): 



f 1 f(x)(^g)(x)dx + j\(x)(3-f)(x)dx = 0, 



(56) 



(57) 



where f E L p (-l,l),p > 1; g e L q (-l,l),q > 1 and — + ~ < 1. We call such a procedure "Tricomi pairing", for which a 
couple of examples are presented in the next proposition. 

Proposition 4.1 (Applications of Tricomi Pairing to Multiple Elliptic Integrals) (a) There are several inte- 
grals that evaluate to the same number |T(^)] 8 /(1287r 2 ): 

J* [k^I-Zs 2 )]^ = y j' 1 DK(*)] 8 dA = 5 plKWfkdk 
= 3 J 1 \K(k)f K[Vl-k 2 ^dk = 2 J^ K ( k ) [^Vl-k^dk = 6j\K(k)] 2 K{yi-k 2 }kdk. 

(b) There are two multiple elliptic integral representations of Apery's constant ((3) = Z^Li n~ 3 : 

((3) = J [2E(k) - K(k)] [K ( y/l - k 2 } ] 3 k dk 

= ^- : |2^ 1 [K(^)] 3 E(v / l-^ 2 )d^+^ 1 [K(v / l-^ 2 )] 3 E(^)(4^-l)d^|. 

Proof (a) If we set fix) = KWd - x)/2)/y/T+x, - 1 < x < 1 and g(x) = 2K(v/(l + x)/2)K( Vd - x)l2), - 1 < x < 1 in Eq.[ 
while recalling the Cauchy principal values given in Eqs.[37]and[52j then we arrive at 



(58) 



K 



1 + x 



K 



1-x 



VT+x J-l \\ 2 ; 



K 



K 



( W 



di 



This instantly rearranges into i / 1 [K(v / l-& 2 )] 3 d& = 3 / 1 [K(£)] 2 K(v / ]T-£ 2 )d& (an identity that Wan conjectured numer- 
ically in (6J] without an analytic proof), which reveals the equivalence between the leading items of the first two lines in 
Eq.gT] 

Writing k = (1 - <f)/d + and employing Landen's transformation K(2\/f/(l + 0) = d + <f)K(0, one has 



/;[K(^] 3 d,=/; 



K A/1 



w 
1+f 



l-f r 1 
— - = 2 / 
l + < Jo 



writing & = (1 - 0/(1 + and employing Landen's transformation 2K((1 - {)/d + 0) = d + OK(\/l-<f 2 ), one has 



Jo Jo I U + wJ l + < 4j 



(1 + 



Od<=-J [k(Vi-& 2 )] dA+-J [K(Tj)] 3 f;df7, 



where the last step is a trivial substitution cf >-» \Jl-r) 2 . The two simultaneous equations displayed above make it pos- 
sible to eliminate any one among the three quantities J 1 [K(vT^I 2 )] 3 d^, $\K{k)fdk, flm.k)fkdk, and determine 
the ratio between the two remaining numbers. Thus, we have verified the chain of identities in the first line of Eq. [57] 
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(cf. [Ref. 6, Eq. 29]). The relations within the second line of Eq.[57]can be likewise established by successive Landen's 
transformations, as detailed in the first paragraph of [Ref. la, p. 139]. 

As we recall from Eq.l26lthat 



the verification is complete, 
(b) InEq.[5lweset 



f(0 = 2K 
g(0 = K 



I /l+fl 


K 








IV 2 J 



-1<<<1; 



+ K 



K 



f 1 + ^ 



E 



{ W 



K<<1, 



to arrive at the following identity: 

j^KtAOK^l-/? 2 ] j - [K(Vl-^ 2 ]] 2 + K(^)E(^) + K(Vl-^ 2 )E(Vl-^ 2 )|^d^ 
= -^jfK^)] 2 - [K(v / l-^ 2 )] 2 |[K(^)K(^l-fe 2 ]+K(v / l-^ 2 )E(fe)-K(fe)E(\/l-^ 2 )]^d^. 
By symmetry under the substitution k >-* Vl-k 2 , the left-hand and right-hand sides of the equation above simplify to 
^ 2K(yfe)K(\/l-fe 2 )|-[K(Vl-^ 2 )] 2 + 2K(Vl-^ 2 )E(Vl-^ 2 )|/fedA! 



and 



- jT* | [K(k)] 2 ~ [K^Vl-k 2 ] ] 2 j [K ( /i^ 2 ) E(« - K(«E ( v 7 !^ 2 ) ] k dk 
■ 2^ [k(^1-^ 2 ]] 2 [K(v / l-^ 2 )E(^)-K(^)E(v / l-^ 2 )]^d^, 



respectively. Upon rearrangement and application of Legendre's relation K(\/t)E(\/ 1 - 1) + K(\/l - t)E(\/i) 
K(Vt)K(VI^t) = n/2 [Ref.[H item 8.122], we obtain 

0= J K(^)[k(v / 1-^ 2 )] &d/e + J [K(Vl-^ 2 )] [K(v / l-^ 2 )E(^)-3K(fe)E(\/l-^ 2 )]^d^ 



= -2 j' 1 K(^)[K(v / l-^ 2 )]%d^ + ^ 1 [k(^1-^ 2 ]] 2 4K(v / l-^ 2 )E(^)-y 



&d&. 



Meanwhile, setting k = (1 - 0/(1 + and applying Landen's transformations, one can verify that 

l3 



^ [K(^l-/j 2 )] 3 E(/z)d/e = 2^ 1 [K( < 0] 3 



Eh/1-^ 2 +fKh/l-< 



so the penultimate equation becomes 



= 2j [K(A0] 3 E(\/l-fe 2 ]dfe + J [k(^1-^ 2 |] E(&)(4/j-l)d/j- —] [k^I-Zz 2 ]] kdk. 

To verify that fQ[K(Vl-k 2 )] 2 kdk = f| /^[P-i^x)] 2 dx = 7<T(3)/4, one may refer to a standard formula (see [Ref.[lS 
item 7.112.3] or [Ref.Hl item 2.17.19.2]) 



[\p v (x)?dx=-?— 
J-i 2v+l 



2sin 2 (v;r) £° 1 
1 t {v+l + n) 2 



, veC\(Z <0 u{-l/2}), 



and take the v — «■ -1/2 limit. This concludes the proof of Eq.[5 
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4.2 Tricomi Transform of P v (x)P v (-x) 
The formula in Eq.[52]can be rewritten as 

- 1 2P_ 1/2 (OP-i/ 2 (-<0 



r 

1 



ji(x - 



df = -{[P_i /2 (x)] 2 - [P_i/ 2 (-x)] 2 }. 



This is not accidental. In the next proposition, we will generalize such a Tricomi transform relation to Legendre 
functions of arbitrary degree v. 



Proposition 4.2 (Tricomi Transform of P v (OPv(-<D) For any v e C \ Z, we have 

J-l 7l{ 



2P v (OPv("0 , - [Py(x)] 2 " [Pv(-X)] 2 

-d<f = - 



r(x - 



sin(v7r) 



For n e Z>o, i/iere is an identity 



r 1 [p^ r 1 iP-n-mf 

J-i2(x-0 J-i 20b- f) 



d^P^QJx). 



(59) 



(60) 



Proof To verify Eq. [59j it would suffice to demonstrate that 

[P v (x)] 2 -[P v (-x)] 2 







sin(v7r) 



J-l 7I(. 



t(x - 



= — — - f P/(x){[P v (x)] 2 -[P v (-x)] 2 }dx+- f Q ? (x)P v (x)P v (-x)dx, WeZ> . 

Sin(V7T)J-l 71 J-l 



(61) 



Here, in the last line, we have used the Parseval identity of Tricomi pairing (Eq.l56D. along with the Neumann integral 
representation for Legendre functions of the second kind (cf. [Ref. 26, Eq. 11.269] or [Ref.[I2, Table 1.12A, Eq. 12A.26]): 



Qel 



x) = &\ -f- 

J-i 2(x 



1 P?(Od£ 




, Vxe(-1,1),WeZ> . 



(62) 



For a non-negative even number £, both addends in the last line of Eq. [61] vanish because the integrands are odd 
functions. We may now settle Eq.[6T]for odd numbers £ by induction. For £ = 1, we use Eq.[16]to compute 



— f 

sin(v7r) J-i 



Pi(x){[P v (x)] 2 -[P v (-x)] 2 }d 



1 r 1 

x — f 

sin(v7r) J-i 



x{[P v (x)] 2 -[P v (-x)] 2 }dx 



4v(v+l)sin(vjr) x— l-o 



lim (l-x / ) — 



dx 



1 , d 

lim (1-x 2 ) — 

4v(v+l)sin(v7r) -i+o + dx 



., 2 d{[P v (x)] 2 -[P v (-x)] 2 } 

(1-x ) 

dx 

„ 2 d{[P v (x)] 2 -[P v (-x)] 2 } 
(1-x ) 



dx 



4sin(v7r) 
v(v + l);r 2 



Meanwhile, we may check that 

4 pi 4 pi 

Q 1 (x)P v (x)P v (-x)dx= - I 

71 J-l 7rJ-l 

= -*Cp v i X )p v i- X )d X - 2 -C log ^ d 

ttJ-i 7rJ-i 
_ 4 r 1 1 
7i J-i 4v(v 



X, 1 + x 

l+-log- 

2 1-x 



P v (x)P v (-x)dx 



4v(v+l) dx 



4v(v+l) dx 

,„ 2^d(Pv(x)P v (-x)) 
(1-X ) 

dx 



(l-x z )— 
dx = 



dx 
4sin(v7r) 



,„ 2 d(P v (x)P v (-x)) 
(1-x )- 



dx 



+ 4v(v + 1)(1 - x 2 )P v (x)P v (-x) di 



v(v + l);r 2 



Thus, Eq.[6T]holds for £= 1. Noting that the Legendre functions of the first and second kinds satisfy similar recursion 
relations, namely, 

(2ju + 1)(1 - x 2 ) ^^ = + l)[P^_i(x) - P M+ i(x)]; + l)xP M (x) = (^ + DP^+itx) + pP^-iix); 

9 dQ,,(x) 

(2ju + 1X1 - x 2 ) ^ = n(fi + l)[Q M _i(x) - Qft+iix)}; (2/x + DxQ^x) = (// + l)Q M+ i(x) + nQ^-iix), 
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we can deduce 



v(v+l)(£ + l) f 1 „ wi ,. . 2 2 



sin(v7r) 
v(v+l)£ 
sin(v7r) 



p e+1 ( X ){[p v (x)r - [p v (- X )Y}dx + 



4v(v+i)(f + i) r 1 



4y(y 4. \\f rl 

P / _ 1 (x){[P v (x)] 2 -[P v (-x)] 2 }dx+ — Q^_i(x)P v (x)P v (-x)dx 

n J-i 



f Q e+1 {.x)P v 



(x)P v (-x)dx 



v ( v + me + 1) r 1 xW{[Pv(x)]2 _ [Pv( _ x)] 2 }dx + 4v(v + iX2/ + i) f xQdx)Pv(x)Pv( _ x)dx , 

J-l 71 J-l 



sin(vjr) 



as well as 



AV ^ + } ] f V + DPe + i(x) + £P e - 1 (x)]{[P v (x)] 2 - [P v (-x)] 2 }dx = 4v(v + 1) f 1 xP/x){[P v (x)] 2 - [P v (-x)] 2 }dx 



= _ f 1 p, (x ) A 1(1 _ x 2) A f (1 _ x2) d{^v(x)] 2 [PA x)] 2 } j + 4v{v + m _ x 2 )[pv(x)] 2 _ 4v(y + _ x 2 )[Pv(x) ?\ dx 

J-i dx { dx [ dx J 



dx 



(l-x z ) 



2 d{[P v (x)] 2 -[P v (-x)] 2 } 



dx 



+ Av{v + 1)(1 - x z )[P v (x)] z - Av{v + 1)(1 - x z )[P v (x)] 



2] dP^(x) 
J dx 



da; 



+ • 



8[l + (-l/]sin 2 (vjr) 



) fa-* 



= £(£ + 1) 

J-l 

8[l + (-l) / ]sin 2 (v7r) 
+ ^ 



V/ ^ dttP^lM^-x)] 2 ) dx + W + l)v(v + l) £ [p ^_ i(x) _ p ^ +i(x)]{[Pv(x)]2 _ [Pv{ _ x)f]dx 



dx 



which leads to a recursion relation 



and 



{£ + l) z [(/ + l) z - (2v + l) z ] P /+1 (x){[P v (x)] z - [P v (-x)] z } dx 



- ^V 2 - (2v+ l) 2 ] f^P e -i(x){[P v (x)f - [P v (-x)] 2 }dx = - 



8[l + (-in(2^+l)sin 2 (v7r) 



4v(v+ 1) 
2^ + 1 



[(£ + DQe+iM + £Q e -i(x)]P v (x)P v (-x)dx = 4v(v + 1) 



J ^xQAx 



)P v (x)P v (-x)dx 



J-i dx ( dx 
Xil^ X ^dx 



„ 2 d(P v (x)P v (-x)) 
(1-x )— 



+ 4v(v + 1)(1 - x z )P v (x)P v (-x) }• dx 
1 dQ^(x) 



dx 

+ 4v(v + 1)(1 - x 2 )P v (x)P v (-x) 



= £(£ + 1) 

2[l + (-l) / ]sin(v7r) 



dx 



2^+1 



j dx 



■dx 



2x<i(Pv(x)Pv(-x)) 

(1-x ) 

dx 



which brings us another recursion relation 

(/+ 1) 2 [(^+ l) 2 -(2v+ l) 2 ] f\ e+ i(x)P v (x)P v (-x)dx 

2[l + (-l) / ](2^ + l)sin(v7r) 



'£< 



-^V 2 -(2v+l) 2 ] / Q/_i(x)P v (x)P v (-x)dx : 



Here, we have used the identities 
lim (l-x 2 )A 



JC-.-1+0+ 



dx 



2(l-x 2 )P v (x) d ^ 
dx 



8 -^pl; P,(l) = (-l/P,(-l) = l,V^Z a0 



to take care of boundary contributions to the integral concerning P((x), and resorted to the limit behavior 

,■ „ 2x2 d QM( x )d(P v (x)Pv(-x)) 2sin(v7r) 

hm (l-x z r — ; ; = 

i-o+ dx dx n 

z.odQJx) d(P v (x)P v (-x)) 2cos(ti7r)sin(v7r) 

hm (1-x) — - = 

sc^-i+o+ dx dx n 
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for the integration by parts involving Qe(x). Therefore, the last line of Eq. [61] satisfies a homogeneous recursion 

(£ + 1) 2 [(£ + if - (2v + 1) 2 ] I — !— [ 1 P (+1 (x){[P v {x)f - [P v (-x)f] dx + - [ 1 Q e+ i(x)P v (x)P v (-x)dx j 

( sin(v7r) J_i it J -i J 



x\ , £eZ> 



= £ 2 [£ 2 - (2v + l) 2 ] I — |— f 1 P^_i(x){[P v (x)] 2 - [P v (-x)] 2 } dx + - f * Q f _i(x)P v (x)P v (-x)d 3 
(sin(vjr)J-i 7rJ_i 

and the truthfulness of Eq.[6T]for £ = 1 entails any scenario with a larger odd number £. This completes the verification 
ofEq.HUfor veC\Z. 

For a fixed x £ (-1,1), both sides of Eq. [59] represent continuous functions of v, so we may handle Eq.[60]by investi- 
gating the v — n limit where n e Z>o. Suppose that n is even and non-negative, then we have 

v [Py(x)] 2 - [P v (-x)] 2 on , _ P v (x)-P v (-x) on , _ cos(v7r)P v (x)-P v (-x) 

hm = 2P„(x) hm = 2P„(x) hm 

v—n sin(v7r) v— n sin(vjr) v— n sin(vnr) 

= -P„(x)limQ v (x) = -P„(x)Q„(x)= -P„(-x)Q„(x). 

7T v—n jr n 

If we start from an odd and positive n instead, we will end up with 

,. [Py(x)] 2 " [P v (-X)] 2 oT1 , N1 . P V (X)+P V (-X) „„ , xl . C0S(V7T)P V (X)-P V (-X) 

hm — — = 2P„(x) hm — — = -2P„(x) hm — — 

n sin(vi) v— n sm(vjr) v— « sin(vsr) 

= - -P„(x) lim Q v (x) = --P„(x)Q„(x) = -P n (-x)Q n (x). 

71 v—n 7r JT 

Hence, we have proved 

2P/i(^)Pn( — £) t i- ^rj, 2P_ tt _i(£)P_ft_i( — £) ( > , * 4 

^ — d<f = ^> - d<r = (-l) ra -P„(x)Q n (x) = -P n (-x)Q n (x), neZ> , 

J-l 7T(x-<) J-l 7T(x-<) 7t 71 

which is equivalent to Eq. [60] In fact, Eq. [60] is not particularly surprising. It is just a special case of the stronger 
statement that (cf. [Ref. 26, Eqs. 11.280 and 11.281] or [Ref. 32, Table 1.12A, Eq. 12A.27]) 



1 Pn(Qp(Q 

2(x-0 



d<f = Q n (x)p(x), degp(x)<rc, (63) 



where p(x) is any polynomial whose degree does not exceed n. 



A key step in the proof of Proposition ^. llta)| hinges on the identity 



/ o 1 [K(>^]^* = aj[W)fK(>^)d*, „ lP dx = s£ p -^g-** dx . 

This result can be further generalized, which in turn gives rise to a family of identities satisfied by multiple elliptic 
integrals involving the product of three complete elliptic integrals (of the first and second kinds). 



Corollary 4.3 (Tricomi Pairing of P(2n+i)/2(xVV 1 + x) For any integer n £ Z, we have 

(2n + l)/2(Q d£ _ Q(2tt + l)/2* 
y/T+l 2(x-0~ VT+X 



9 r 1 p^mio dt = q^^x) Vx 1} (64) 

J-l X 1 " 



and 



11 [P(2 n + 1)/2(X)] 3 , „ f 1 P(2 ra + l)/2(x)[P( 2n+ l)/2(-^)] 2 



J-l vT+x J-l 



vTTx J-i \/i + x 

particular, for n-0, one has 



dx. (65) 



j' 1 [2E(Vl-^ 2 )-K(Vl-^ 2 )] 3 d^ = 3^ 1 [2E(^)-K(^)] 2 [2E(Vl-^ 2 )-K(Vl-^ 2 )]d^. 
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Proof We might recall from the definition of the Legendre functions of the second kind (Eq. [3]l that 

Ql2»+lVa<*) = (-D n+1 ^P(2n + l)/2(-x). 

Therefore, the two formulae in Eqs. 1371 and 1381 correspond to Eq. [64] for the cases where n — -1 and n — 0. We now 
prove E q . l64l inductively for positive integers n e Z>o. Suppose that Eq.|64]is true in the range neZn[-l,m] for some 
non-negative integer m e Z>o, then we may rely on the recursion relations 

2m + 3 2m + 1 

2(771 + l)xP( 2m + l)/2(x) = P(2m+3)/2W+ £ P(2m-l)l2W, 

2m + 3 2m +1 

2{m + l)acQ(2OT+l)/2(«) = „ Q(2m+3)/2(*)+ = Q(2m-1)/2W, 



for the computation 

2m + 3^ r 1 P( 2m+3 )/2(0 d< 



U Ji 



1 P( 2m+ i)/2(0 fd£ 2m + l T 1 P (2m -i)/2(0 d<f 



= 2( " , + 1 ^/-i ^/irfVx-O 2 38 j_! X TT7 2(.v-c) 



J-i n/T 



„, Q(2m + l)/2(^) , 
; 2(m + 1 Yr (to + 1) 

Vl + x 

2m + 3 Q(2m+3)/2(x) 



^ P (2m+1)/2(Q d c 2m +1 Q(2m-l)/2(x) 

s/T+l 2 vT+x 

1 P(2m+1)/2(«D , , 2m + 3 Q(2m+3)/2(^) 



T 1 P(2m+1 
J-l 0^ 
f 1 P(2m + 



(m + l) ^z±^Zl d ^.. . 

2 v^+x J-l ^/T+t 2 v^+x 

To justify the last step, we use the Mehler-Dirichlet integral representation for Legendre functions (Eq. Q} to compute 



(m+l) f P (2m +i)/2(cos0) 
Jo 



sin0d0 
VTTcosB 



-(m+l) 



(m+l) 



r \ 2 f 

JO 7t Jo 
JO 71 Jfl 



cos((m + 1)/3) 
o ^(cos/S- cos 8) 



dp 



sin0d0 



x/TTcos© 



sin0d0 



V^cos/S-cosS) \/l + COS0 

V2(m + 1) / cos((m + l)/3)d/3 = 0, Vm £ Z, 
Jo 



cos((m + l)/3)d/3 



where the interchange of the order of integrations is permitted by absolute integrability and the Fubini theorem. This 
eventually shows that the truthfulness of Eq.[64]for n e Z n [-l,m] entails the corresponding scenario where n- m+l: 



1 P( 2to+3)/2(£) _ Q(2m+3)/2(*) 

^1+7 2(x-0~ y/T+X ' 



thus Eq.[64]holds for all positive integers n e Z>o. The conclusion extends to negative integers by symmetry: P(2n-l)/2 : 

P(-2n-l)/2 and Q(2n-l)/2 = Q(-2n-l)/2- 

Before showing Eq.[65l we recall from Eq.l59lthat 



f 1 1 1 Nn 2 - P (2n + l)/2(OP(2n+l)/2(-f) , r rD , .,2 I D , ,,2 

J (-IF cx d< = [P(2„ + l)/2te)] " [P(2«+l)/2(-^)] 



-1 < X < 1. 



1 < X < 1 



7T(x - 

This formula can be combined with an equivalent formulation of Eq. [64] 

p f 1 P(2»+i)/2(Q d<f = n+1 P(2n+i)/2(-*) 

+ x 

in an application of Tricomi pairing (Eq.[56). The result in Eq.|65]is then immediate. 

We follow up with some additional examples involving the product of four elliptic integrals in the integrands. 
Corollary 4.4 (Another Application of Tricomi Pairing) We have an integral identity 

sin(2v7r)cos(v7r) 



J x[P v (x)] 3 P v (-x)dx = 



(2v+l) 2 n 



veC\{-l/2}, 



(66) 
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which leads to 



71 

2 



9^ 



l/2(x)] P_i/2(-x)dx 



—j (l-20[K( v ^)] 3 K(v / l r 7)dt=--7 / (l-20[K(vT : 0] 3 K(v / Odf, 

7T* Jo 7T 4 Jo 



(67) 



_i/ 3 (x)] P_i/ 3 (-x)dx 



216 f 1 (l-p 2 )p(2 + p) 



6_ p d-p- 

T 4 Jo 1 



y/En* Jo 1 + 2p 

72 ri(l-p a )p(2+p) 



1-2 



27p 2 (l + p) 2 



7 2 r 1 
In 4 Jo 



2\/2 

71 



y/Zn 4 Jo 1 + 2p 

i/4(x)] 3 P_i/4(-x)dx 



1-2 



4(l + p + p 2 ) 3 
27p 2 (l + p) 2 



K 



I lpH2 + p)\ 



l + 2p j 



K 



p 3 (2 + p)' 
l + 2p 



dp 



4(l + p + p 2 ) 3 



K 



p 3 (2 + p)' 
l + 2p ; 



K 



( pH2 + p)\ 



l + 2p 



dp, 



(68) 



£> 

_ 32\/2 r 1 l-2u 
7i 4 Jo (1 + A) 2 

64^ f 1 l-2w 



K 



2^ 
1+ \/w 



K 



1 - \/w 



k V 1 + 



dw 



K 



( j i-VU 

l + Vu 
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>4V2 r l 

l^Jo (1+v^) 2 

J x[P_i/ 6 (x)] 3 P_i/ 6 (-x)dx 

_ 54 T 1 f(l - 0(1 + 1)(2 - 0(1 - 20 
~^Jo 

54 r 

?r 4 Jo 



K 



2V 7 " 

i+ V" 



du, 



(69) 



(l-i + * 2 ) 3 
54 r 1 t(l- 0(1 + 0(2 -0(1 -20 
o 



(l-* + i 2 )' 



[K(v / F)] 3 K(vT^7)df 
[K(vT^7)] 3 K( v / Odi. 



(70) 



Proof For integer degrees v e Z, the left-hand side of Eq. [66] represents the integration of an odd function over the 
interval [-1, 1], hence vanishing. This is consistent with the corresponding behavior on the right-hand side. 

We now turn our attention to the scenarios where veC\(Zu(- 1/2}). From the Tricomi transform formula in Eq. |591 
one may readily deduce the following relation 



1 2fP v «)P v (-f) 



t(x-<) J-i 



J-l IF 



1 2[x - (x - OFvlW-O , , [Pv(x)] 2 - [Pv(-x)] 2 



7T(x - £) 



-d£ = x- 



sin(vjr) 



n J-i 



(OP v (-0df. 



Here, according to Eq. I19|9m „)1 we have 



«)P v (-Od<: 



2cos(v7r) 
2v+l ' 



so we may combine the formula 



1 2^P v (0Pv(-0 



J-l 7r(x 



d<f = x 



[P v (x)] 2 - [P v (-x)] 2 4cos(v;Q 
(2v+l)jr 



- sin(v7r) 
with Eq.[59]for the implementation of the following Tricomi pairing: 

** 2P v (0P v (-{)sin(vjr) 



j ^ xP v (x)P v (-x){[P v (x)] 2 - [P v (-x)] 2 } Ax = j ^ xP v (x)P v (-x) \& j ^ 

= - J 1 2fPv(f)Pv(-S)sin(v70 

J-i 



7l{x - 



dx 



ti(x - 



d£ 



)Pv(-x){[P v (x)] 2 - [Pv("X)] 2 }dx + 



P v (x)P v (-x)dx 
2sin(2v7r) 



(2v+l)jT 



J P v (x)P v (-x)dx. 



After rearrangement, we obtain 



/l />i 2sin(2v7r) 
x[P v (x)] 3 P v (-x)dx = 2 / xP v (x)P v (-x){[P v (x)] 2 - [P v (-x)] 2 }dx = — — - \ P v (x)P v (-x)dx, 
-i J-i (2v + l)7r J-i 
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which entails the claimed identity in Eq. [66] 

The left-hand side of Eq. [66] extends to be a continuous function in v e C. Taking the v - 
Eg. [671 The special cases v = -1/3,-1/4,-1/6 correspond to Egs. l68l[69l and[70l 

4.3 Tricomi Transform of (1 + x) y - n P v (x),n £ Z> ,Re(v-rc) > -1 



- 1/2 limit, one arrives at 



In Eq. [64] we evaluated the Tricomi transform of P(2 ra +i)/2(x)/\/l + x,n £ Z, which extends the Neumann integral 
representation of Legendre functions Qe tor £ £ Z>o (Eq.[62l. 

In this subsection, we present a generalization of Eq. [64] to Legendre functions of arbitrary complex- valued degree 
v e C. As we shall soon see, the integral representation of P v in Eq.[T]is virtually not used in the evaluation of a related 
Tricomi transform. Our proof in the next proposition is of algebraic/combinatorial (as opposed to geometric/analytic) 
flavor, which hinges on recursion relations rather than manipulations of multiple integrals. Apart from the Parseval- 
type identity for Tricomi transforms (Eq.l56ll. we shall also make use of the Hardy-Poincare-Bertrand formula (see [Ref. 
ml §4.3, Eq. 4] or [Ref.lH, Eq. 11.52]): 



ST Vf{STg) + g(tf-f)] = (3-f)(3-g) - fg, 

which applies to the scenarios where f e IP(-1, l),p >l;ge LH-1, l),q > 1 and | + | < 1. 

Proposition 4.5 (Generalized Neumann Integrals) Denote by {-v} the fractional part of -v e (-oo, 1) \ ; 
have the generalized Neumann integral representation for Legendre functions of the second kind: 



(71) 



then we 



p v (0 as 



Qv(x) 



-i(l + { - v) 2(x-0 (l + x) ( - v> 
Furthermore, for n e Z>o and Re(v- n) > —1, one has 



-1 < x < 1. 



(72) 



l 2(x-<) 



-■(l + x) v ~ n Q v (x), -1<x<1. 



(73) 



Proof The following is a variation on the proof of Eq.(64] Recalling the recursion relations for Legendre functions 

(2p + l)xP^(x) = (^ + DP^+iix) + iiP^-xix); {2\x + DxQ^x) = (ju + DQ^+xix) + ^Q^_i(x), 

we may deduce the following identities for {-p} = {-v}: 



= (2p+l) 



Q^ix) /-i p„ +1 (0 as 



a+x) { - v] 

xQ^ix) 



= (2p + l)x 



(l + x) { - vl 
Q»{x) 
(l + x) { - v] 



( l + f )(-v} 2(x-0 

P^O (S-x + x)dt 
( l + £)(-v> 2(x-0 
P^O d( 



Qfi-i(x) 
(l + x) { " v > 



1 Pu-l(0 d{ 



{ ~ v) 2(x-0 



(1 + {_V} 2(x-0 



2\i+\ 



i: 



(i + H 



Here, the last integral can be computed with the aid of a standard integral formula [Ref. ll3l item 7.127]: 

Recr > -1, 



I 



Si \fj t) / \ i 2 a+1 \Y{o+l)f 
(1 + x) P M (x)dx = 



r(o-+ju+2)r(i + (r-ju) 



which entails the vanishing identity 



I 



and the homogeneous recursion 

Qn+i(x) 



(2jU+l)x 



(l + x) ( - vl 
Q M (x) 



i(l + { - v 

i p^iO at 



d£ = 0, (j + (-v)eZ >0 



(74) 



(75) 



(l + x) { - vl 



(l + { - v] 2(x-0 
P»iO d$ 



(l + f)f-v) 2(x-0 



(l + x) { " v > 
ju + {-v}eZ>o. 



i P^_!(0 d£ 



2(x-0 



(76) 
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Thanks to the recursion in Eq. [76] it would suffice to verify Eq. [72] in two distinguished scenarios: -1 < v < and 
0< v< 1. 

For -1 < Rev < 0, we consider the expression 

p v (x):=-(l + x) v Q v (x)+?-&> (\l + O v Pv(0 rx , -1<*<1, (77) 
2 J-x 7i(x-0 

and investigate its recursive properties in the next three paragraphs. 

First, by definition of p v (x) in Eq.l77land an application of Eq.[74j we have 

P - v -x(x) = -(l + x)- v - 1 Q- v -x(x)+ f 1 (l + 0" v " 1 f-v-i(0 7^ =r , (78) 

2 J-x n(x-0 

(1 + x)p v (x) - ™ V + o! ]2 = " (1 + ^) V+1 Qv(^) + f ) v+1 P v (0-*L-. (79) 

T(2v + 2) 2 J_i jr(x-<f) 

Taking the last pair of equations (Eqs. [78] and [79} as inputs for the Hardy-Poincare-Bertrand formula (Eq. [TP , we 
obtain 



f 1 P v (OWv(0 + Q-v-i(0]-^-r + ^ [ 1 ^v(Ot(l + 0" v Pv(0 + (1 + O v+ Vv-i(0] 

2 J-i 7T(x-4) 2 J-i 



7T(X - <f) 



2 

= Q v (x)Q_v_i(x) - — [P v (x)f + (1 + x)" v Q_v-i(a:)pv(a:) + d + x) v+1 Q v (x)p- v -x(x) + (1 + x)p v (x)p_ v _i(x). 

According to the relation between Q v and P v (Eq.[3) and the Tricomi transform of P v (0P v (-£) (Eq.[59l, the last equation 
reduces to 



({)[(! + <)" v Pv(0 + (1 + f) v+ p-v-i(f)] 



7T(X - 



= (1 + jc) v Q_ v _i(x)p v (x) + (1 + x) v+1 Q v (x)p_v_i(x) + (1 + x)p v (x)p- v -x(x). (80) 
As we set v - + irj,ri eU in the equation above, we obtain 

l^/Vi^ORe^ (81) 

a functional relation that will become useful later. 

Next, in view of the shifted moment formula for Tricomi transforms (cf. [Ref.[2|| Eq. 11.217]): 

» f 1 (1+ f = a + f - - La + *) m f 1 (i + sr-^nodt, e z> , (82) 

J-X 7l(X-i) J-XTl(X-i) TT^O J-l 

we may extend Eq.l79lto 

(i+x) n+1 p v {x)- \ £d+*r [\i+o v+n - m p v (Odt 
= -a+x) v+n+1 Q v (x)+^&> (\i+o v+n+1 P v (0^-r v M n) ) 

2 J-x n(x-0 
Pairing up Eqs. [78] and [Hl^] in the Parseval identity (Eo.l56t. we obtain 

- n n ] : f 1 (l + ,, Pv«)Pv(-Odf = J f 1 (l + n P v (O[Qv«) + Q-v-l(O]d^ 
2sm(v7r) J-i 2 J-i 

= - J f 1 (i+o"p v (o[(i+o- v p v (o+(i+0 v+1 p-v-i(0]d^+^ £ / >1 (i+xr v - 1+m p v (x)dx f d+er+^PyCOdf. 

2 J-i 4 J^qJ-x J-x 

We assert that the equation above entails 

J\l + O n P v (0[d + 0~ v p v (S) + (1 + O v+1 p- v -x(Om = 0, -1< Re v < 0,n £ Z> (83) 
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owing to the presence of an identity of convolution type for -1 < Rev < 0,re £ Z>o: 

'2 



- ^ : f 1 (l + f) B Pv(?)Pv(-Odf=7 £ f 1 (l + ^r v - 1+m P v (x)d a; f 1 (l + f) V+n " m i , v«)d<. 
2sin(v7r)J-i 4^ J_i J-i 



(84) 



It takes some effort to prove such a combinatorial relation as Eq. [84] As we consider the left-hand side of Eq. |84j we 
have 



2sim 



^ 2 r 1 „ 

in(vu) J-i 



(OP v (-f)d{=- 



2sin( 



n 2 r 1 
in(v7r) J-i 



(l + OPv(OPv(-fldf = 



7T 2 COt(V7T) 

2v+l 



according to the formula Jj^PvCOPvC-Odf = 2cos(v;r)/(2v + 1) established in Eq. I19|9m ^1 and the Appell-Legendre 
differential equation (cf. Eq.[I"6l> 



d£ 



d£ 



2 d(P v (QP v (-Q) 
d<f 



+ 4v(v + 1)(1 - r )P v (f)Pv(-0 f + 4v(v + DfP v (OP v (-0 = 



allows us to deduce the following recursion for re > 2 via simple integration by parts: 
M n := f\l + O n P v (OP v (-Od( 



-L 
-L 



(l + O ra_1 Pv(f)Pv(-Odf- 



(1 + <) ™-i Pv(<)Pv( _^ )d ^ + . 



4v(v 
re-1 



— r 

£ 



(1 + B_1 T7 ^(l-D-TT 



d£ 



d<T 



s d<f 



+ 4v(v + ixi - n \ p v (OPv(-odi 



+ (»-!) j" 1 



4v(v+ 1) 
(1 + 0"" z d - ^ z )Pv(OPv(-^)d^ 



(1 + 0"" 2 (e 2 « 2 - 2<fre 2 + re 2 + ( 2 n + 4fre - 5re + 4)P v «)Pv(-<)d< 



(re-1) 3 (2re-l)[2v(v+l)-re(re-l)l (re- l)[-4v(v+ l) + n(n+ 1)] „ 
M„_ 2 + — M„_! + — M n . 



v(v+l) 2v(v+l) 
In other words, the sequence 



4v(v+ 1) 



2 sin' 



n 2 r 1 
in(v7r) J-i 



(l + f)"P v (OPv(-Odf = 



2sin(v7r) 



Mn 



is recursively characterized by 

7T 2 COt(V7r) 

re(n - 2v - l)(re + 2v + 1)L„ + 2(2re - l)[2v(v + 1) - re(re - l)]L„_i + 4(re - 1) 3 L„_ 2 = 0, re £ Z> 2 . 
Meanwhile, the right-hand side of Eq. [84] can be recast in terms of gamma functions by virtue of Eq.[74l 

2" +1 [T(m - v)] 2 [T(-m + n + v+ l)] 2 



(85) 



fV + xr v - 1+m Pv(*)dxf\l + v+ "- m Pv(0d£=7 £ 



3, T(m + l)r(m - 2v)r(-m + re + DlX-rei + re + 2v + 2) 



Direct computation reveals that 



Po=Pi 



7T COt(V7T) 

2v+l 



We then show that {P n |re £ Z>o} satisfies the same three-term recursion as the sequence {L n \n £ Z>o}. Clearly, the 
generating function 



n , s V r, » 71 V 2 f - v [r(^-v)]^x^ 
G R (x):= > P„x =-•<> 



2 I co om+v+1 



[T(rez + v+l)r x 



2 „m 



7T 2 COt(V7T) 

2v + l 



2-PL 



-V, -V 

-2v 



2x 2 Pl 



m=0 T(m + 2v + 2) m\ 

v+ l,v+ 1 
2(v+l) 



2x 
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is the product of two hypergeometric functions. It can be verified that the following second-order ordinary differential 
equation 



4(2x - l)x 

has two linearly independent solutions 
f 1 (x) = X - &v+1)/2 2 F 1 



2 d /Cx) .,. „ &f(x) m ^ 



da: 2 



■ + 4(4x- l)x- 



dx 



■ + [2x + (2v+in/Xx)=0 



-v,-v 
-2v 



2x] and /-2(x) = x (2v+1)/2 2Fi 



V+ l,v+ 1 
2(v+l) 



2x 



According to a special case of Appell's theory all the three products [/i(x)] 2 , /i(x)/2(x) and [^(x)] 2 satisfy a third- 
order differential equation 



(2x - lYx 



2 3 d 6 g(x) 2 A S(x) r _ 2 o .^dgd) , , . :i 



dx 3 



■ + 3(2x-l)(4x-l)x z 



dx 2 



■ + 4x[7x z + (2v^ + 2v - 3)x - v(v + 1)] - 



dx 



+ 4x(x + v + v)g(x) = 0. 



In particular, the generating function Gr(x), being a constant multiple of /i(x)/2(x), must solve the differential equation 
above. This fact can be paraphrased as 

x— [x— - 2v - 1) [x— + 2v+ 1) G R (x) + 2 (2x— - 1 



dx \ dx 



dx 



V dx 



2v(v+ l)-x — x 1 

dx \ dx 



(xG R (x)) + 4\x—-l\ (x'G r (x)) = 0, 



which is equivalent to the recursion relation n(n-2v-l)(n+2v+l)R n +2(2n-l)[2v(v+l)-n(n-l)]R n -i+4(n-l) 3 R n -2 = 
0,n e Z>2. As the sequences {L n \n e Z>o) and {R n \n e Z>o) on both sides of Eq. [84] share the same recursive characteri- 
zation (Eq.[55]l, they must be identical. Consequently, the truth of Eq. [Ml implies that of Eq.|8"31 
After this, with the sequence of relations 



J P. 1+^(0(1 + O n Re [(1 + f)H«J p _ i^^)] df = 0, re e Z> ,t7 1 



in hand, and the fact that P_i , > for -1 < <f < 1, we may deduce 

2 T "/ 

Re[(l + O»" i, P_j +jl ,(O]=0, -l<f<l 
by continuity. Going back to Eq.[81j we then have 



(86) 



= Re (l + x)s-»iQ _ 4 _ i -Cx)p i +i _C*) + (l + x^p i +i <x) 



Im [(1 + x)s- i "p_ i +ifJ (*)] {-|P_ i +i;? (x)tanh(r/7r) + Im [(1 + x)5" ir 'p_ i + } , -1< x < 1. 



2" -2+"? 

Combining the results from Eqs.l8"6"landl87l we see that 



(1 + x) v p v (x) |--P v (x)cot(v7r) + (l + x) v p v (x)]=0, -1<x<1 



(87) 



(88) 



holds for Rev = -1/2. By analytic continuation, the vanishing identity in Eq.[88]remains true for -1 < Rev < 0. As the 
two factors on the left-hand side of Eq. [88] are both continuous in x e (-1,1) and analytic with respect to v in the strip 
-1 < Rev < 0, only one of them can vanish identically. From the Neumann integral representation of Qo (Eq.l62D. we 
know that lim v ^o Pv(x) = 0, -1 < x < 1, so the factor 



-— P v (x)cot(v7r) + (l + x) v p v (x) 

cannot be identically zero. By the dichotomy in Eq.[88l we must have p v (x) = 0, thus reaching the conclusion that Eq.[72l 
holds for -1 < v < 0. 

To settle Eq.[72]for < v < 1, we note that the expression p v (x) defined in Eq. 1771 extends to be an analytic function 
of v for Rev > -1. By the principle of analytic continuation, we have the following identity for Rev > -1: 



(l + x) v Q v (x): 



7i r 1 
2 J-i 



Jl(x - O 



1 < X < 1. 
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Meanwhile, we may enlist Eq.[75]for the computation of 



2 J-i 7r(x-0 2 J-i n(x-0 2 J-i 



which consequently confirms Eq.[72]for < v < 1. 

Now that we have established the Tricomi transform formula in Eq. [72] we may also extend the arguments in the 
last paragraph to prove Eq.[73]by induction. ■ 

Remark One may well recognize that Eq.[73]incorporates the classical result in Eq.[63]as a special case: 

- 1 P n (Op(0 



•L 



1 2(x-0 



■dZ = Q n (x)p(x), degp(x)<neZ> . 



While the proposition above yields formulae that extend the Neumann integral (Eq. [62l . there was not a new proof 
for the latter. In fact, in our proof, we have presented the generalized Neumann integral in Eq. [72] effectively as 
a "functional square root" (in the Hardy-Poincare-Bertrand sense) of the Tricomi transform formula for P v (<f)Pv(-<0 
(Eq.l59D. the latter of which was proved in Proposition [4j2] using the original Neumann integral (Eq. I62l>. 
One may find the following formula, tabulated without proof or reference, as item 2.17.3.5 in 12811 : 

" P iX + 0) P v {~) dx = -(a + y) v Q v [-), Rev > -l;y « [-a,aT 
J - a x-y \a> \a> 

If we conveniently ignore a (probably misprinted) factor of 2, we may consider the n — scenario of Eq. [73] as an 
"analytic continuation" of the equation above. However, unfortunately, the definition of Q v (z),z e C \ [-1,1] has not 
been disambiguated in |28], and the precise meaning of "y t [-a, a]" is not clear. As a result, it is fairly difficult to verify 
the numerical/analytic reliability of item 2.17.3.5 in (280, m particular for complex- valued v. Nonetheless, we are under 
the impression that such a tabulated formula is reparable, and perhaps there exists a purely analytic proof of Eq. [73l 
without using recursions or a detour to the Tricomi transform formula for P v (f )P v (-£) (Eq. I59l>. 

Had one been given Eq. [73] as a known input verifiable by purely analytic means, the proof of the generalized 
Neumann integral (Eq. [72l would be made as concise as that of Eq. [64] Consequently, one would also obtain the 
Tricomi transform formula for P v (f)P v (-f) (Eq. I59l> as a corollary, thus dispensing with the recursive proof devised for 
the latter in Proposition [472] In a certain sense, Eqs. [59] and [72] are two equivalent formulations of the same fact. 

Unlike the geometrically-motivated proof of Eq. [64] or the classical analysis behind the Neumann integral formula 
(Eq.l62|l. the demonstration of the generalized Neumann integral formulae in this section apparently draws very little 
on the analytic properties of the fractional degree Legendre functions. However, we point out that there are indeed a 
couple of analytic ingredients vital to the generalized Neumann integral formulae, as briefly explained below: 

(1) The integral formula in Eq.[74l 

This is undeniably a consequence of the integral representation for P v (Eq.Q}. We have the following computations 
for Re A < 1: 



2 prccosf cOS^Hli^ 

v/2(cos0-O 
l r n 
(A)f(§-A) Jo '< 



r 1 P^«)d< r 1 2 f™f cos ^ 2 r 71 r C( 

J-i2(l + A ~ J-i nlo J2(cosB-n d ^ (T+O 1 ~ n Jo J-i 



V2(cos^-0( 1 + 



(2/i + l)j3 Ja 
cos d p 



cos — cos dp, 



0rsin(Ajr)r(A)r(|-A)Jo (l + cos^) A 2 

where the interchange of the order of integrations is justified by absolute integrability and the integration in ( e 
(-1, cosp") can be completed by hypergeometric techniques. For A e (-oo, 1/2) u (1/2, 1) and arbitrary /u, one can use the 
knowledge of hypergeometric functions to evaluate 



I 



cos — cos dp = 



o (l + cos^) A 2 2 r sin(2A7r)r(2A-l)r(-A-^+l)r(-A + ^ + 2) 

as well. In other words, a self-contained proof of the generalized Neumann integral formulae would require the integral 
representation for Legendre functions (Eq.Q} at its initial stage. 

(2) The applicability of "Tricomi pairing" 

To justify various applications of the Parseval identity (Eq. l56l l and the Hardy-Poincare-Bertrand formula (Eq. I71l > 
during the "algebraic/combinatorial proof" of the generalized Neumann integral formulae, we need an a priori estimate 
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for the integrability property of the function p v (x) as defined in Eq. [77] One may set such a bound leniently, say, by 
using the triangle inequality and the boundedness of the operator 3~ : L p (-1, 1) — ► L p (-1, 1) (1 < p < +oo) to put down 



JjPvCxFdx) " <[j\l+x) v Q v (x)\ p dx^ P + | (1^1^(1 + O v Pv(0- 



(x-0 

r 1 \ 1/p n ( r 1 ^ Vp 

J ^\a + x) v Q v (x)\ p &xj +C P -[J |(l + x) v P v (x)| p dx 



P \Vp 
dx 



where the power p is chosen to ensure convergence of both integrals in the last line, and the finite constant C p is 
selected per the operator norm of 3~ : L p (-1, 1) — ► L p (-1, 1). □ 

We wrap up with some applications of the generalized Neumann integral (Eq. 172} to multiple elliptic integrals 
whose integrands contain six factors of Legendre functions. 

Corollary 4.6 (Tricomi Pairing) (a) For veC\(Zu {-1/2}), we have 

3/ P v (-x)[P v (x)]°dx=[2sin z (v7r) + 3] / [P v (x)P v (-x)] 3 dx, (89) 



(90) 



J' 1 i P v (-x)[P v (x)] 5 dx= [2sin 2 (v7r) + 3] fiPvi 
3 f 1 x 2 P v (-x)[P v (x)] 5 dx=[2sin 2 (v7r) + 3] x 2 [P v (x)P v (-x)] 3 dx + 16C ° S (v7r)sm {vn \ 

J-l J-l (2V+1) 3 7T 2 

(b) The following identities hold: 

/l g fi 5\/3 

i [P_ 1/6 (x)P_ 1/6 (-x)] 3 dx = - J i P_i/ 6 (-x)[P_ 1/6 (x)] 5 dx = -y- J i [P_ 1/6 (-x)] 2 [P_ 1/6 (x)] 4 dx, (91) 

6 j' i X 2 P_i/6(-x)[P_ 1 / 6 (x)] 5 dx-15v / 3j' i X 2 [P_i/ 6 (-x)] 2 [P_i/ 6 (x)] 4 dx + 20 J x 2 [P_ 1/6 (x)P_i/ 6 (-x)] 3 dx = 0. (92) 

Proof (a) In the Hardy-Poincare-Bertrand formula (Eq.\TQ, we take f(() = g(() = 2P v (<f)Pv(-0, and invoke Eq.H to 
deduce 



J-l tt(x-<) sin(v7r) { sm(vn) 



Now, applying the Parseval identity (Eq.[56jl to the pair of Tricomi transforms in Eqs.[59]and[93] we obtain 

6 f 1 P v (x)P v (-x) | [Pvfa)] . ~ [Pvi ~ x)] I dx = 8[ 1 [P v (x)P v (-x)fdx, 
J-i { sm(vn) J J-i 

which rearranges into Eq. [89] 

A similar "Tricomi pairing" procedure for /"(£) = g(0 = 2(P v (£)P v (-0 bring us 



4£ 2 P v (OPv(-O[Pv(<0]-[Pv(-O] 2 Jr 2 |[Pv(x)] 2 -[Pv(-x)] 2 l n m/ln( - 16C0S 2 (V7T) 

7 ^ —< \ d t= X 1 —, \ \ ~ 4X IP V (X)P V (-X)] +— T-n-n, 

7r(x-<) sin(v7r) ( sin(v7r) J (2v+lrjr z 



so a subsequent application of the Parseval identity (Eq.[56]l results in 

[P v (x)] 2 - [P v (-x)] 2 ] 2 , o f 1 2r „ , _ , „ 3 l 64cos 3 (v;r) 

(2v+1) 3 tt 2 ' 



f 1 9 ( [Py(x)] 2 - [P v (-x)] 2 1 f 1 , o 

6/ x 2 P v (x)P v (-xH . , , — — ^dx = 8l x 2 [P v (x)P v (-x)] 3 dx- 

J-i ( sin(v7r) J J-i 



which implies Eq. [90] 

(b) For - 1/2 < v < 0, one may directly apply the Parseval identity (Eq. l56l and the Hardy-Poincare-Bertrand formula 
(Eq.[7T]l to two copies of Eq.[72l and deduce 



(l + O 2v Pv(OQ v (Odt = 0, (95) 



(i + 2v P v (OQ v (0-^ = a+xf v \ [Q v {x)f - —[Pyix)f \ , am 
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respectively. A combination of the two equations above leads to another formula: 

(l + x-x + <f)d<f 



J-i x-t, J -I 



x-i, 



= i\+x) 2v+1 \[Q v {x)f-—{P v (x)f\- \ (l + £) 2v P v (OQ,(Odc 



-/;< 



= (1 + x) 2v+1 j[Q v (x)] 2 - — [P v (x)f J . (97) 
Specializing Eqs. [96] and [97] to the case where v = -1/6, we obtain a pair of equations 

.'/> I (l+0~ 1/3 P-i/6(OQ-i/6(0— r ^a + x)- 1/3 \[Q-ye(x)f-^[P- V6 (x)f\, (98) 



X - £ 

,^ . r „ , 

x-£ I 4 



3» / (l + 2/3 P-i/6(OQ-i/6(0-^7 = (l + x) 2/3 nQ_i/ 6 (x)] z --—[P- ye (x)r\, (99) 



which further iterate into 



25» ^(l + ^-^P-i/eCOQ-i/eCojtQ-i/eC^)] 2 - ^[P-i/eGO] 2 } 



x-£ 



,2 

: (1 + x)" 2/3 \ [Q_i/ 6 (x)] 2 - ^-[P_ 1/6 (x)] 2 \ - (1 + x)- 2/3 7T 2 [P_ 1/6 (x)Q_ 1 / 6 (x)] 2 , (100) 



23» / (1 + <) 4/3 P_ 1/6 (OQ-i/6(0 <! [Q-i/e(0] z - ^[P-i/eCO] 



4 ] x-$ 

i 2 > 2 

= (1 + x) 4/3 1 [Q_i/ 6 (x)] 2 - ^[P_ 1/6 (x)] 2 L - (1 + x) 4/3 7r 2 [P_ 1/6 (x)Q_ 1/6 (x)] 2 . (101) 

Pairing up Eqs.l99"land ll00l in the Parseval identity (Eq.[56l>. we obtain a vanishing identity 

2 J P-P + (6Pt-15v / 3P 3 P-+20P 2 P 2 .)d < f = J P_P + (6Pl - lSv^P+P 3 + 40P+P 2 . - lSv^P+P- + 6P^)d^ 0, (102) 

where we have used the short-hand notations P+ = P_i/6(±0- Meanwhile, from a special case of Eq. [89] we have 
Gj\P-P 5 + d£ = 7/_ 1 1 P 3 _P 3 d£, so the last identity also implies 5^ f^P-Pi d£ = g/^P^P 3 d<f. This proves Eq.[91] 

If we pair up Eqs.l98land ll01l in the Parseval identity (Eg. [561. we obtain another vanishing identity 

j 1 P_P + (6P 4 - lSv^P+P 3 + 40P 2 P 2 - lSv^P+P- + 6P 4 )(1 + f)df = 0, 

which is not substantially more informative than Eq. 11021 This is because we have 

J P-P + (QP 4 _ - 15VIP+P 3 + 40P 2 P 2 . - lSv^P+P- + 6P\)$ d£ = 

by symmetry. Bearing in mind the equation above, we may consider applying the Parseval identity (Eq. [56]l to Eqs. l99l 
and llOll and simplify the result 

J 1 P-P + (6Pt - 15\/3P + P 3 _ + 40P 2 P 2 . - lSv^P+P- + 6P 4 )(1 + O 2 d( = 

into 

J P-P + (6Pt ~ ISv^P+P- + 2QP 2 + P 2 _)f d<f = 0, 
which proves Eq.[92] ■ 
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5 Discussion and Outlook 



In our evaluations of some generalized Clebsch-Gordan integrals (Eqs. I26I29I I. the expressions involving Euler's 
gamma function actually correspond to certain special values of complete elliptic integrals as well. For example, 
the knowledge of K(l/v^) = [T(\)fli^\/n) [Ref. HE item 8.129.1] allows us to recast the formula [T(±)] 8 /(1287r 2 ) = 
f^KiVl-k 2 )] 3 dk into the following form: 



r a -i ] 4 =r[f 

Jo vT^sVl-Cs 2 ^) J o Jo 



dt 



lO VT^Wy/l-(l-k 2 )t 2 



dk. (103) 



Here, both sides are integrations of algebraic functions over algebraic domains, which qualify them as members in the 
"ring of periods" denned by Kontsevich and Zagier [1]. It is generally believed that identities for periods can be proved 
by "algebraic means" Oil, namely, relying on nothing else than additivity of the integral, algebraic change of variables, 
and the Newton-Leibniz-Stokes formula. However, the analytic proof we produced for the generalized Clebsch-Gordan 
integrals does not fall into such a category: we have invoked Bessel functions, which are "exponential periods" fl]. We 
would like to see purely algebraic evaluations of the generalized Clebsch-Gordan integrals when the choice of degree 
v leads to special values of complete elliptic integrals. In particular, it could be interesting to search for potential 
connections to high degree modular equations and the Chowla-Selberg theory. (After communicating the first version 
of this manuscript to Prof. Jonathan M. Borwein in January 2013, I was sent a preview of a forthcoming book 13311 
that contains a sketched proof for Eq. 11031 using lattice sums and modular forms, totally independent of the methods 
presented in my work. Later on, James G. Wan also wrote me about his plan to give a fuller account for the proofs of 
his own conjectures in a joint work with Rogers and Zucker, currently available as [34].) 

The spherical methods in this work (Legendre functions and spherical rotations) enable us to handle a large variety 
of multiple elliptic integrals, but they are by no means a cure-all. The methods of Bessel moments 10, 51 Eh, geometric 
transformations of Watson type 13], hypergeometric summations (gl still play fundamental roles in our quantitative 
understandings for integrals over elliptic integrals. 

With a synthesis of various techniques, it is sometimes possible to compute integrals over the product of more than 
three Legendre functions of the same degree v e C, such as 



. .. 4 , .. 2 sin 4 (7rz)[y< 2 >(z + 1) + i// 2 >(-z) + ggftgjj , 2) , d 3 

,(x)rdx=hm — —5-^ , where y^'(z) := -—;logr(z). (104) 

«— v (2z + lrn i dz d 



One may wish to compare Eq. 11041 to the integrals f^ 1 x[P v (x)] 3 P v (-x)dx,v e C (Eq. [66]! that reduce to elementary 
functions. Clearly, special cases of Eq. 11041 bring us some interesting evaluations of multiple elliptic integrals. For 
example, we have the following integral representations for Apery's constant: 



n C n C 2n r 

C(3) = -— x[P- VZ (x)fdx = -—\ x[P- Vi (x)fdx = -—\ x[P- 1/6 (x)fdx 
z4d J-i lboj-i lesyj-i 

as well as a critical scenario involving ((5) = J_^ =1 n~ 5 : 

f(5) = -^r C x[P-y 2 (x)fdx = J- C{2t-l)[K{Vt)fdt. 
372 J-i 93 Jo 



References 

[1] Maxim Kontsevich and Don Zagier. Periods. In Bjorn Enquist and Wilfried Schmid, editors, Mathematics Unlim- 
ited — 2001 and Beyond, pages 771-808. Springer, Berlin, Germany, 2001. 

[2] D. H. Bailey, J. M. Borwein, D. J. Broadhurst, and M. L. Glasser. Elliptic integral evaluations of Bessel moments 
and applications. J. Phys. A, 41:5203-5231, 2008. 

[3] I. J. Zucker. 70+ years of the Watson integrals. J. Stat. Phys., 145:591-612, 2011. 

[4] Jonathan M. Borwein, Dirk Nuyens, Armin Straub, and James Wan. Some arithmetic properties of random walk 

integrals. Ramanujan J., 26:109-132, 2011. 
[5] Jonathan M. Borwein, Armin Straub, James Wan, and Wadim Zudilin. Densities of short uniform random walks. 

Canad. J. Math., 64:961-990, 2012. (with an appendix by Don Zagier). 
[6] James G. Wan. Moments of products of elliptic integrals. Adv. Appl. Math., 48:121-141, 2012. 
[7] Joseph Kleiber. On a class of functions derivable from the complete elliptic integrals, and connected with Legen- 

dre's functions. Messenger Math., 22:1-44, 1893. 



38 



[8] E. W. Hobson. The Theory of Spherical and Ellipsoidal Harmonics. Cambridge University Press, Cambridge, UK, 
1931. 

[9] Srinivasa Ramanujan. Modular equations and approximations to n. Quart. J. Math., 45:350-372, 1914. 
[10] Bruce C. Berndt. Ramanujan's Notebooks (Part V). Springer- Verlag, New York, NY, 1998. 
[11] Bruce C. Berndt. Ramanujan's Notebooks (Part III). Springer- Verlag, New York, NY, 1991. 

[12] Robert Fricke. Die elliptischen Funktionen und ihre Anwendungen. Erster Teil: Die funktionentheoretischen und 

analytischen Grundlagen. B. G. Teubner, Leipzig, Germany, 1916. 
[13] I. S. Gradshteyn and I. M. Ryzhik. Table of Integrals, Series, and Products. Academic Press, Burlington, MA, 

7th edition, 2007. (translated from the Russian by Technica Scripta, Inc., edited by Alan Jeffrey and Daniel 

Zwillinger). 

[14] Wilhelm Magnus, Fritz Oberhettinger, and Raj Pal Soni. Formulas and Theorems for the Special Functions of 
Mathematical Physics, volume 52 of Grundlehren der mathematischen Wissenschaften. Springer, New York, NY, 
3rd enlarged edition, 1966. 

[15] G. N. Watson. A Treatise on the Theory ofBessel Functions. Cambridge Mathematical Library. Cambridge Univer- 
sity Press, Cambridge, UK, 2nd edition, 1995. 

[16] M. L. Glasser. Definite integrals of the complete elliptic integral K. J. Res. Natl. Bur. Stand. B: Math. Sci., 
80B:313-323, 1976. 

[17] P. Appell. Memoire sur les equations differentielles lineaires. Annates scientifiques de I'E. N. S., 10:391-424, 1881. 
[18] Harry Bateman. Higher Transcendental Functions, volume I of Bateman Manuscript Project. McGraw-Hill, New 

York, NY, 1953. (compiled by Arthur Erdelyi, Wilhelm Magnus, Fritz Oberhettinger, Francesco G. Tricomi, David 

Bertin, W. B. Fulks, A. R. Harvey, D. L. Thomsen, Jr., Maria A. Weber and E. L. Whitney). 
[19] Elias M. Stein and Rami Shakarchi. Complex Analysis, volume II of Princeton Lectures in Analysis. Princeton 

University Press, Princeton, NJ, 2003. 
[20] F. W. J. Olver. Asymptotics and Special Functions. Computer Science and Applied Mathematics. Academic Press, 

New York, NY, 1974. 

[21] D. S. Jones. Asymptotics of the hypergeometric function. Math. Meth. Appl. Sci., 24:369-389, 2001. 

[22] A. P. Prudnikov, Yu. A. Brychkov, and O. I. Marichev. Integrals and Series, volume 2: Special Functions. Gordon 

and Breach Science Publishers, New York, NY, 1986. = A. n. IIpynHHKOB, K). A. BptraKOB, O. PI. MaptraeB. 

HnmezpaAu u psrdbi. Tom 2. Cnev,uajibHue §ynKv,uu. Hayica, MocKBa, CCCP, 1983 (Translated from the 

Russian by N. M. Queen). 

[23] A. Gervois and H. Navelet. Some integrals involving three Bessel functions when their arguments satisfy the 
triangle inequalities. J. Math. Phys., 25:3350-3356, 1984. 

[24] E. Beltrami. Intorno ad un teorema di Abel e ad alcune sue applicazioni. Reale Istituto Lombardo di Scienze e 
Lettere, Rendiconti, 13:327-337, 1880. 

[25] Paul F. Byrd and Morris D. Friedman. Handbook of Elliptic Integrals for Engineers and Scientists, volume 67 of 
Grundlehren der mathematischen Wissenschaften. Springer, Berlin, Germany, 2nd edition, 1971. 

[26] Frederick W. King. Hilbert Transforms (Volume 1), volume 124 of Encyclopedia of Mathematics and Its Applica- 
tions. Cambridge University Press, Cambridge, UK, 2009. 

[27] Alexander Apelblat. Table of Definite and Infinite Integrals, volume 13 of Physical Sciences Data. Elsevier, Ams- 
terdam, The Netherlands, 1983. 

[28] A. P. Prudnikov, Yu. A. Brychkov, and O. I. Marichev. Integrals and Series, volume 3: More Special Functions. 
Gordon and Breach Science Publishers, New York, NY, 1990. = A. II. npy/mnKOB, IO. A. Bpmmkob, O. M. 
Mapn^CB. HnmezpaAu u pfrdu. Tom 3. Cnev,uajibHue tf>yHKv,uu. JIonojmumejibHbie ZAaeu. HayKa, MocKBa, 
CCCP, 1986 (Translated from the Russian by G. G. Gould). 

[29] W. N. Bailey. A reducible case of the fourth type of Appell's hypergeometric functions of two variables. Quart. J. 
Math., 4:305-308, 1933. 

[30] Francesco Giacomo Tricomi. Integral Equations. Dover Publications, New York, NY, 1985. 

[31] Elias M. Stein and Guido Weiss. Introduction to Fourier Analysis on Euclidean Spaces, volume 32 of Princeton 
Mathematical Series. Princeton University Press, Princeton, NJ, 1971. 

[32] Frederick W. King. Hilbert Transforms (Volume 2), volume 125 of Encyclopedia of Mathematics and Its Applica- 
tions. Cambridge University Press, Cambridge, UK, 2009. 

[33] J. M. Borwein, M. L. Glasser, R. C. McPhedran, J. G. Wan, and I. J. Zucker. Lattice Sums: Then and Now, volume 
150 of Encyclopedia of Mathematics and Its Applications. Cambridge University Press, August 2013. 

[34] M. Rogers, J. G. Wan, and I. J. Zucker. Moments of elliptic integrals and critical L-values. arXiv: 1303.2259vl 
[math. NT], 2013. 



39 



